Chapter 


Failure Theories 


5.1 Preliminary Discussions 


The modes of failure discussed in Chapter 2 were presented as an overview of the possible 
mechanisms which trigger a failure. Assessing failure requires an accurate and thorough 
stress analysis in which the largest stresses (principal stresses) that occur at a critical point 
in the component under evaluation must be estimated. In addition to defining the principal 
stresses, one must also be aware of the type of material they are working with. Ductile and 
brittle materials will typically fail in different manners, and their failure is predicted by 
different theories. Similarly, when fatigue is being considered, the overall strength of the 
material is reduced as the number of load cycles increases. 


5.2 Multiaxial States of Stress and Strain 


Under the most complicated loading conditions, a machine part may be subjected to forces 
and moments produced by various combinations of direct axial loads, bending loads, direct 
or transverse shearing loads, torsional loads, and/or surface contact loads. As noted in 4.2, 
no matter how complicated the geometry of the part or how many different external forces 
and moments are being applied, the force system can always be resolved into just three re¬ 
sultant forces and three resultant moments, each defined with respect to an arbitrarily cho¬ 
sen x-y-z coordinate system. The origin of the selected x-y-z coordinate system may be 
placed at any desired critical point within the part. The most complicated state of stress that 
can be produced in a small elemental volume of material at the critical point is the triaxial 
state of stress 1 ', any triaxial state of stress may be fully specified relative to the chosen 
x-y-z coordinate system by specifying the three normal stress components a x , a y , a v and the 
three shearing stress components r xy , r yz , and r 7X . 

Principal Stresses 

It can be shown 2 that if the selected x-y-z coordinate axes, together with the elemental vol- 
ume of material, are rotated in three-dimensional space, about the fixed origin, a unique 
new orientation may be found for which the shearing stress components vanish on all faces 
of the newly oriented elemental volume. This unique orientation, for which the shearing 
stress components on all faces of the elemental cube are zero, is called the principal ori¬ 
entation. For the principal orientation, the three mutually perpendicular planes of zero 


] See Figure 4.1. 2 For example, see ref. 1, Ch. 4. 
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shear are called principal planes , and the normal stresses on these principal planes (planes 
of zero shear) are called principal normal stresses , or just principal stresses. 

The importance of principal stresses lies in the fact that among them will be the 
largest normal stress that can occur on any plane passing through the point, for the given 
loading. The principal stresses are usually designated a h a 2 , and a 3 . There are always 
three principal stresses, but some of them may be zero. 

There is a second special orientation of coordinate axes that may be found, for which 
the shearing stresses on the faces of the rotated volume element reach extreme values. Of 
the three extreme values, one is a maximum, one is a minimum, and the third is a minimax 
value. This orientation is called the principal shearing orientation , and the shearing 
stresses on these three mutually perpendicular principal shearing planes are called princi¬ 
pal shearing stresses. The principal shearing planes may have nonzero normal stresses act¬ 
ing on them, depending upon the type of loading, but normal stresses on planes of 
principal shear are not principal stresses. 

The importance of principal shearing stresses lies in the fact that among them will be 
the largest or maximum shearing stress that can occur on any plane passing through the 
point, for the given loading. The principal shearing stresses are usually designated r h r 2 , 
and t 3 . 

Stress Cubic Equation 

Depending upon whether its material behaves in a brittle or a ductile manner, failure at the 
governing critical point of a machine part is dependent upon the principle normal stresses, 
the principal shearing stresses, or some combination of these. In any event, it is important 
for a designer to be able to calculate principal normal stresses and principal shearing stresses 
for any combination of applied loads. To do this, the general stress cubic equation* may be 
employed to find the principal stresses cr 1 , <j 2 , and a 3 as a function of the readily calculable 
components of stress a x , a y , a v r xy , r yz , and r zx relative to any selected x-y-z coordinate sys¬ 
tem. The general stress cubic equation, developed from equilibrium concepts, is 

cr 3 - a 2 (a x + a y + a z ) + a{a x a y + a y a z + <t z ct x - t% - t 2 z - t%2) 


°v4>) = 0 (5-1) 


(tT x dy (T ? “I - '2T x yTy-,T. 


xy T yz T zx cr x T yz ^y^zx 


Since all normal and shearing stress components are real numbers, all three roots of 
(5-1) are real; these three roots are the principal normal stresses oq, a 2 , and a 3 . 

It is also possible to find the directions of principal stress vectors (and principal shear¬ 
ing stress vectors) if necessary. 4 

Furthermore, it can be shown that the magnitudes of the principal shearing stresses 
may be calculated from 



(T 2 (J 3 

2 


(5-2) 



°~3 (J \ 

2 


(5-3) 



(J 1 (J 2 
2 


(5-4) 


3 For example, see ref. 1, p. 97. 4 See, for example, ref. 1, 
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To summarize, if loads and geometry are known for a machine part, a designer may 
identify the critical point, arbitrarily select a convenient x-y-z coordinate system, and cal¬ 
culate the resultant six stress components cr x , cr y , a v t^, r yz , and r^. Equations (5-1) 
through (5-4) may then be solved to find the principal normal stresses and the principal 
shearing stresses. This procedure works for all cases; if the case is biaxial, one of the prin¬ 
cipal stresses will be zero, and if it is uniaxial, two principal stresses will be zero. 


Example 5.1 Multiaxial States of Stress; Principal Stresses 

A hollow cylindrical 4340 steel member has a 1.0-inch outside diameter, wall thickness of 
0.25 inch, and length of 30.0 inches. As shown in Figure E5.1A, the member is simply sup¬ 
ported at the ends, and symmetrically loaded at the one-third points by 1000-lb loads. The 
tubular bar is simultaneously subjected to an axial force of 5000-lb tension and a torsional 
moment of 3000 in-lb. For the critical point at midspan, determine (a) the principal stresses 
and (b) the principal shearing stresses. 

Solution 

a. Since the vertical transverse loads and end supports are symmetric, the vertical com¬ 
ponents (z-components) produce pure bending, with the maximum bending moment 
extending over the entire middle span. Therefore bending stress (tensile) will be max¬ 
imum at the lower extreme fiber all along the center span between the 1000-lb trans¬ 
verse loads. The axial tensile force produces uniform axial (x-component) tensile 
stress throughout the wall. The torsional shearing stress produced by torque T will be 
maximum at the outer fibers throughout the tube. Based on these observations, the 
midspan critical point A is depicted in Figure E5.1B. The resultant state of stress at 
critical point A due to bending, tension, and torsional shear is shown in Figure E5.1C. 

The stress components a x and T xy may be evaluated as follows: 

Me F a 

(T , =-I-— 


and 


T 


xy 


Ta 

J 


From the dimensions of the tube 


~ df) 
64 


tt(1.0 4 - 0.5 4 ) 
64 


0.046 in 4 



Figure E5.1A 

Hollow bar loaded by forces and torques. 


Figure E5.1B 

Critical point location. 


Figure £5.1C 

State of stress at critical 
point A. 
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Example 5.1 
Continues 


Therefore 


[<1°00>< 10 >](t) 5000 

0.046 + 0.589 


117,200 psi 


and 



0.092 


16,300 psi 


For the state of stress at critical point A depicted in Figure E5.1C, the stress 
cubic (5-1) reduces to 

o- 3 - a 2 a x + a-(-Tly) = 0 


giving three roots 



These may be rewritten as 



03 = 0 


Substituting for cr x and r xy , the principal stresses become 


& l = 


117,200 

2 


117,200 V 


(16,300) 2 


= 58,600 + 60,825 = 119,425 psi 
cr 2 = 58,600 - 60,825 = -2225 psi 

cr 3 — 0 


b. From (5-2), (5-3), and (5-4) and using results from above the principal shearing stresses 
become 


k 2 1 
k 3 | 


-2225 - 0 


0 - 119,425 


= 1113 psi 


* 59,713 psi 


119,425 - (-2225) 


= 60,825 psi 
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Example 5.2 Principal Stress 

It is desired to assemble a cast-iron hub to a steel shaft using an interference fit for 
retention. The system is subjected to a vertical load P and a torque T as shown in 
Figure E5.2 (a). At the critical location the stress field is given by the stresses shown in 
Figure E5.2 (b) 




Figure E 5.2 

Cast-iron hub fit to a steel shaft and associated state of stress at the critical location. 


Solution 

The stress cubic equation is given as 

ct 3 — a 2 1 1 + dll ~ h = 0 


where 


I\ — &X &Z 

i-2 ^y^z ^x&z Try ^yz T zx 

It, — (T x (TyCr z + '2T X yTy 7 T zx CT x Ty Z (JyT zx (J Z T X y 


are the stress invariants. The above values can be written in matrix form as 



<r x 

T tv 

Txz 


20 

7.5 

15.5 

(T = 

Txy 


T yz 

= 

7.5 

15 

-12.5 


_ r xz 

T yz 

<*z_ 


15.5 

-12.5 

-30 


Using any number of programs, such as Maple, Matlab, etc., the above can be solved for 
the roots of the polynomial equation or the matrix equation can be solved for the eigen¬ 
values, which are the principal stresses of the system. These are given as 


15 MPa 


o-i = 26.22 MPa, 


o - 2 = 16.94 MPa, 


0-3 = -38.16 MPa, 
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Mohr's Circle Analogy for Stress 

For biaxial states of stress the stress cubic equation (5-1) degenerates to the form 

a 3 - a 2 (a x + a y ) + cr{a x a y - r 2 y ) = 0 (5-5) 

if stress components lie only in the xy plane. If stress components lie only in the yz plane 
or the zx plane, the stress cubic yields the same formal expression with solely yz subscripts 
or solely zx subscripts, respectively. The principal stress solutions for (5-5) are (see 
Example 5.1) 




<r 2 


a 3 


a x + a y 
2 

(T X + (Ty 
2 

= 0 


+ 



(5-6) 

(5-7) 

(5-8) 


Examination of (5-6) and (5-7) suggests that these equations are formally the same as 
the equation of a circle plotted on the cr-r plane, since the equation of such a circle on an 
xy plane has the form 


(x ~ h) 2 + (y - k) 2 = R 2 


(5-9) 


where h and k are the x-y coordinates of the center of the circle, and R is the radius. An anal¬ 
ogy (refer to the discussion of the “membrane analogy” in 4.4) was developed by Mohr 5 in 
1882 in which he successfully postulated that a circle plotted on the <t-t plane could be 
used to represent any biaxial state of stress because the governing equations are formally 
the same for a circle as for a biaxial state of stress. Such plots have come to be known as 
Mohr’s circles for stress. For example, a biaxial state of stress in the xy plane (see principal 
stress solutions in (5-6), (5-7), and (5-8)) may be represented by a circle 6 with center on the 
(j-axis of a cr-r plot if in (5-9) x and y are replaced by a and r respectively, and 


h 


d x J r a 


2 


k = 0 


R 



(5-10) 

(5-11) 

(5-12) 


Figure 5.1 shows such a plot. To successfuly utilize Mohr’s circle as a stress analysis tool 
the following sign conventions should be adopted: 

1. Normal stresses should be plotted as positive if tensile, and negative if compressive. 

2. Shearing stresses should be plotted as positive if they produce a clockwise (CW) cou¬ 
ple, and negative if they produce a counterclockwise (CCW) couple. (This convention 
is used only for Mohr’s circle applications.) 

The plotting procedure, illustrated in Figure 5.1, is to plot point A (tensile a y and CW 
Tjcy couple), then point B (tensile a x and CCW couple), then pass a circle through A and 
B. The location of center C and magnitude of radius R may be established geometrically 


5 See ref. 1. 6 For example, see ref. 1, or ref. 8, Ch. 10. 
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-t 





(i a ) Biaxial element at 
critical point. 


(CW) 



Figure 5.7 

Mohr’s circle for a biaxial 
state of stress in the x-y 
plane. 


( b ) Mohr's circle for xy biaxial state of stress shown 
in (a); note that cr 3 = 0 for this biaxial case. 


as shown, and principal stresses and a 2 pinpointed where the circle intersects the cr-axis 
(shear is 0 on the cr-axis). The principal shearing stress t 3 may also be found as shown (see 
(5-4)). The third principal stress, <r 3 = 0 , should not be ignored; plot a 3 at the origin of the 
cr-r plane and then construct two additional Mohr’s circles, one through <r 2 - cr 3 and an¬ 
other through - a 3 . From these circles r, and t 2 may be found as shown. If the maxi¬ 
mum shearing stress for this state of stress is needed for failure prevention calculations, it 
will be the largest of r 1? r 2 , and t 3 ; for the case illustrated in Figure 5.1 this is r max = r 2 . 
Note specifically that r max is not r 3 , which is based on the initially plotted solid circle. If 
cr 3 = 0 had not been plotted, the dashed circles would have been missed, and, therefore, 
the true maximum shearing stress would have been missed. 

It is also possible to find the orientation of the principal axes from the Mohr’s circle 
construction, using the expression 


2r xy 

26 = tan 1 - (5-13) 

o-y ~ (t x 

where the double angle 20 is shown in Figure 5.1(b), measured/rc>m the ray CA to the prin¬ 
cipal stress ray a x . This corresponds to the rotation angle 6 in Figure 5.1(a) measured in 
the same direction as 20, from the y-axis to the principal axis 1. A consequence of this pro¬ 
cedure is that principal axis 2 will be perpendicular to axis 1 in the x-y plane, and princi¬ 
pal axis 3 will be mutually perpendicular to axes 1 and 2. 


Example 5.3 Mohr's Circle for Stress; Principal Stresses 

Using the Mohr’s circle analogy, solve the problem posed in Example 5.1 to find 

a. The principal stresses 

b. The principal shearing stresses 
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Example 5.3 
Continues 


Z, 


X 



y 


y 



T xy = 16,300 psi 


c j x = 117,200 psi 


Figure E5.3A Figure E5.3B 

State of stress at critical point A. Biaxial xy element from Figure E5.3A. 


Figure E5.3C 

Mohr’s circle constructed for biaxial 
element shown in Figure E5.3B. 



Solution 

Reviewing the analysis of Example 5.1, the state of stress at the critical point A is shown 
in Figure E5.3A, repeated from Figure E5.3C. 

Also, the magnitudes of cr x , a y , and r xy were determined in Example 5.1 to be 

(T x = 117,200 psi 
cr y = 0 

r xy = 16,300 psi 

a. Using the Mohr’s circle analogy to find principal stresses, the xy plane of Figure 
E5.3A may be redrawn as shown in Figure E5.3B. Mohr’s circle may next be con¬ 
structed as shown in Figure E5.3C by plotting points A and B as shown, passing a cir¬ 
cle through them and finding principal stresses and <x 2 semigraphically. The third 
principal stress, cr 3 = 0, is plotted at the origin. 

From Figure E5.3C, the following values may be found graphically 

_ _ 117,200 

a x = C + R = ---+ 60,825 = 119,425 psi 

_ _ 117,200 

a 2 = C - R = --- 60,825 = -2225 psi 


and for this biaxial state of stress 


cj 3 = 0 
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These values are the same as those found in Example 5.1 by solving the stress 
cubic equation. 

b. To find the principal shearing stresses, the remaining two (dashed) Mohr’s circles are 
constructed, and the principal shearing stresses are found to be, from the radii of the 
three circles, 


|Til = ^2-3 = —=1113 psi 
119 425 

k 2 l =*,- 3 = -^—= 59,713 psi 
|r 3 1 = R\- 2 = 60,825 psi 

These values are the same as those shown in Example 5.1. From these values it may 
be noted that 

Tmax = t 3 = 60,825 psi 


Strain Cubic Equation and Principal Strains 

Just as stress is an important measure of loading severity, strain may also be an impor¬ 
tant parameter in assessing failure potential in a machine part. Strain is the term used 
to define the intensity and direction of the deformation at a given critical point, with re¬ 
spect to a specified plane or set of planes passing through the critical point. Analogous 
to state of stress, the state of strain at a point may be fully defined by six components 
of strain relative to any selected x-y-z coordinate system: three normal strain compo¬ 
nents and three shearing strain components. The normal strain components are usually 
designated s x , e^, and s z , and the shearing strain components y xy , y yz , and y^. Just as the 
state of stress may be completely defined in terms of three principal stresses and their 
directions, the state of strain may be completely defined in terms of three principal 
strains and their directions. The principal strains may be determined from a general 
strain cubic equation analogous to the general stress cubic equation (5-1). The strain 
cubic equation is 


e 3 - s 2 [s x + By + B z ] + 8 


s x s y s y s z S z S x ^ (Try Ifyz Tzx) 

,2 ' 


A (yxyYyzyzx e ;yYzx s zYxy) 


= 0 


(5-14) 


This equation is formally the same as the stress cubic equation (5-1) except that the nor¬ 
mal strains, and the shearing strains divided by two, replace the normal stresses and 
shearing stresses. The three solutions of (5-14) are the three principal normal strains e 1? 

e 2 , and s 3 . 

Mohr's Circle Analogy for Strain 

In some circumstances, a designer may be interested in finding stresses or strains at par¬ 
ticular points in an existing machine part. Such an interest may stem from the need to im¬ 
prove failure resistance, to verify a calculation procedure, or some other requirement. In 
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these cases it is common to bond a strain gage or a strain rosette to the surface at one or 
more critical points. 7 Using the strain gage data to evaluate biaxial solutions of the strain 
cubic equation (5-14), the principal strains and their directions may be found. Use of a 
Mohr’s strain circle to determine the principal strains is common. 

Mohr’s strain circles may be constructed using the techniques just described for 
Mohr’s stress circles, except that normal and shearing strains are plotted on the cartesian 
axes e and y/2 rather than a and r. To utilize the resulting Mohr’s strain circle solutions 
for principal strains, with the objective of estimating stresses at the critical point, relation¬ 
ships between stress and strain must be known. 


Elastic Stress-Strain Relationships (Hooke's Law) 

For a broad class of engineering materials it has been experimentally well established that 
a linear relationship exists between stress and strain so long as the material is not loaded 
beyond its elastic range. These linear elastic relationships were first presented in the 
seventeenth century by Robert Hooke, and have come to be known as the Hooke ’s Law 
relationships. They are 


and 


e x = ~[<Tx - v(cr y + crj] (5-15) 

e y = ~[(T y - v(a z + O-J] (5-16) 

S Z = - v (<Tx + 0>)] (5-17) 




(5-18) 

(5-19) 

(5-20) 


where the elastic constants are 

E = Young’s modulus of elasticity 
v = Poisson’s ratio 
G = shear modulus of elasticity 

Numerical values of these elastic constants are given in Table 3.9 for several engi¬ 
neering materials. 

7 A strain gage is a device designed to measure the displacement between two points some known distance 
apart. They may be mechanical, optical, electroresistive, piezoresistive, capacitive, inductive, or acoustical in 
nature. The electrical-resistance strain gage is often employed on the free surface of a machine part or speci¬ 
men to find stresses there. It is usual to use an arrangement of three gages to measure three strains at known 
angles to each other, from which principal surface strains and directions may be found. These three-gage arrays 
are called rosettes. To convert strains to stresses, Young’s modulus of elasticity and Poisson’s ratio must also 
be known. 


Stress Concentration 215 


It should be noted that the Hooke’s Law relationships are valid for any arbitrarily 
selected x-y-z cartesian coordinate system orientation, including the principal 1-2-3 orien¬ 
tation. The Hooke’s Law relationships may also be inverted to give 


<7x = S 


[(1 - v)e x + v(s y + 8 Z )] 

(5-21) 

(1 

- v - 2v 2 ) 



a y = m 

E 

[(1 - v)s y + v(e z + ej] 

(5-22) 

(1 

- v - 2v 2 ) 




E 

o..2\ 

[(1 - v)e z + v(e x + gj,)] 

(5-23) 


(1 - v - 2v 2 ) 


5.3 Stress Concentration 


Failures in machines and structures virtually always initiate at sites of local stress concen¬ 
tration caused by geometrical or microstructural discontinuities. These stress concentra¬ 
tions, or stress raisers , often lead to local stresses many times higher than the nominal net 
section stress that would be calculated without considering stress concentration effects. An 
intuitive appreciation of the stress concentration associated with a geometrical discontinu¬ 
ity may be developed by thinking in terms of force flow through a member as it is subjected 
to external loads (see 4.3). The sketches of Figure 5.2 illustrate the concept. In Figure 
5.2(a) the rectangular flat plate of width w and thickness t is fixed at its lower edge and 
subjected to a total vertical force F uniformly distributed along the upper edge. The dashed 
lines each represent a fixed quantum of force, and the local spacing between lines is there¬ 
fore an indication of the local force intensity, or stress. In Figure 5.2(a) the lines are uni¬ 
formly spaced throughout the plate, and the stress a is uniform and calculable as 



In the sketch of Figure 5.2(b) a flat rectangular plate of the same thickness is subjected 
to the same total force F, but the plate has been made wider and notched back to provide 
the same net section width w at the site of the notch. The lines of force flow may be visu¬ 
alized in very much the same way that streamlines would be visualized in the steady flow 
of a fluid through an obstructed channel with the same cross-sectional shape as the plate 
cross section. It may be noted in Figure 5.2(b) that no force can be supported across the 
notch, and therefore the lines of force flow must pass around the root of the notch, just as 
a flowing fluid must pass around the obstruction. In so doing, force-flow lines crowd 


tttttttt 


. 


F 



Figure 5.2 

Intuitive concept of stress concentra¬ 
tion. (a) Without stress concentration, 
(b) With stress concentration. 


(a) 


(b) 
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Figure 5.3 

Some common examples of stress concentration, (a) Gear teeth, (b) Shaft key way. (c) Bolt threads, (d) Shaft 
shoulder, (e) Riveted or bolted joint, (f) Welded joint. 


together locally near the root of the notch, producing a higher force intensity, or stress, at 
the notch root. Even though the net section nominal stress is still properly calculated by 
(5-24), the actual local stress at the root of the notch may be many times higher than the 
calculated nominal stress. 

Some common examples of stress concentration are illustrated in Figure 5.3. 
Discontinuities at the roots of gear teeth, at the corners of key ways in shafting, at the roots 
of screw threads, at the fillets of shaft shoulders, around rivet holes and bolt holes, and in 
the neighborhood of welded joints all constitute stress raisers that usually must be consid¬ 
ered by a designer. The seriousness of the stress concentration depends on the type of load¬ 
ing, the type of material, and the size and shape of the discontinuity. 


Stress Concentration Effects 

For the purposes of studying stress concentration effects, stress raisers may be classified 
as being either highly local or widely distributed. Highly local stress raisers are those for 
which the volume of material containing the concentration of stress is a negligibly small 
portion of the overall volume of the stressed member. Widely distributed stress raisers are 
those for which the volume of material containing the concentration of stress is a signifi¬ 
cant portion of the overall volume of the stressed member. For the case of a highly local 
stress concentration, the overall size and shape of the stressed part would not be signifi¬ 
cantly changed by yielding in the region of the stress concentration. For a widely distrib¬ 
uted stress concentration, the overall size and shape of the stressed part would be subject 
to significant changes by virtue of yielding in the region of stress concentration. Small 
holes and fillets would usually be regarded as highly local stress concentrations; curved 
beams, curved hooks, and eye-and-clevis joints would usually be categorized as cases of 
widely distributed stress concentration. 

With the foregoing definitions made, stress concentration effects may be classified as 
shown in Table 5.1. From the table it may be noted that the effects of stress concentration 


TABLE 5.1 Effects of Stress Concentration 


Type of Stress 
Concentration 

Type of 
Foading 

Type of 
Material 

Type of Failure 

Stress Concentration 
Factor 

Widely distributed 

Static 

Ductile 

Widely distributed yielding 

K t (modified) 

Widely distributed 

Static 

Brittle 

Brittle fracture 

K, 

Widely distributed 

Cyclic 

Any 

Fatigue failure 

K f 

Highly local 

Static 

Ductile 

No failure (redistribution occurs) 

1 

Highly local 

Static 

Brittle 

Brittle fracture 

K, 

Highly local 

Cyclic 

Any 

Fatigue failure 

K f 
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must be considered for all combinations of geometry, loading, and material with the pos¬ 
sible exception of one, namely, the case of a highly local stress concentration in a ductile 
material subjected to static loading. In this case the local yielding is usually negligible and 
a stress concentration factor of unity may often be used. All other cases must be analyzed 
for potential failure because of the effects of stress concentration. The final column in 
Table 5.1 lists K t or K f as stress concentration factors. The factor K t is the theoretical elas¬ 
tic stress concentration factor, which is defined to be the ratio of the actual maximum lo¬ 
cal stress in the region of the discontinuity to the nominal net section stress calculated by 
simple theory as if the discontinuity exerted no stress concentration effect; that is, 

actual maximum stress cr t 

K t = -—--= — (5.25) 

nominal stress cr nom 

It should be noted that the value of K t is valid only for stress levels within the elastic range, 
and it must be suitably modified if stresses are in the plastic range. 

The factor K f is the fatigue stress concentration factor, which is defined to be the ratio 
of the effective fatigue stress that actually exists at the root of the discontinuity to the nom¬ 
inal fatigue stress calculated as if the notch has no stress concentration effect. K f may also 
be defined as the ratio of the fatigue strength of an unnotched member to the fatigue 
strength of the same member with a notch. Thus, the fatigue stress concentration factor 
may be defined as 


effective fatigue stress S N (unnotched) 
f nominal fatigue stress $n( notched) 

Stress concentration factors are determined in a variety of different ways, including direct 
measurement of strain, utilization of photoelastic techniques, application of the principles 
of the theory of elasticity, and finite element analysis. Numerical values for a wide variety 
of geometries and types of loading are presented in reference 9. A few of the more com¬ 
mon cases drawn from the literature are reproduced in Figures 5.4 through 5.12. Data for 
the important case of screw threads are sparse, but reference 9 reports a range of values for 
K t from about 2.7 to 6.7 at the thread root of standard threads. Most (but not all) threaded 
fastener failures tend to occur in the thread at the nut face. 


Multiple Notches 

Sometimes it will be found that one stress raiser is superimposed upon another, such as 
a notch within a notch or a notch in a fillet. Although accurate calculation of the overall 
stress concentration factor is difficult for such combinations, reasonable estimates can 
be made. 8 Figure 5.13(a) illustrates a large notch with a smaller notch at its root. To es¬ 
timate the combined influence of these notches, the stress concentration factor K tX for 
the large notch is determined as if the small notch does not exist. This permits an esti¬ 
mate of the stress a' n at the root of the large notch by multiplying the nominal stress a n 
by K tl to give 


<r' n = K tl a n (5-27) 

Now, assuming that the stress a' n occurs throughout the zone within the dashed line near 
the notch root of Figure 5.13(a), cr' n becomes the nominal stress as far as the small notch 
is concerned (because the entire small notch lies within the a' n field). Next, determining 


’See ref. 6. 
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Figure 5A 

Stress concentration factors for a shaft 
with a fillet subjected to (a) bending, 

(b) axial load, or (c) torsion. (From ref. 2; 
adapted with permission of John Wiley 
& Sons, Inc.) 
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Figure 5.5 

Stress concentration factors for a shaft 
with a groove subjected to (a) bending, 
(b) axial load, or (c) torsion. (From 
ref. 2; adapted with permission from 
John Wiley & Sons, Inc.) 


(a) 



(b) 



(c) 
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rid r/d 

(b) ( c ) 


Figure 5.6 

Stress concentration factors for (a) a shaft with a radial hole subjected to axial load, bend¬ 
ing, or torsion (from ref. 3 with permission of the McGraw-Hill Companies), (b) a shaft 
with a straight parallel key way subjected to torsion (from ref. 4, with permission of the 
McGraw-Hill Companies), and (c) a shaft with an eight-tooth spline subjected to torsion 
(from ref. 2, with permission from John Wiley & Sons, Inc.). 
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(a) 


(b) 


Figure 5.7 

Stress concentration factors for a flat bar with a shoulder fillet subjected to (a) bending or 
(b) axial load. (From ref. 2; adapted with permission from John Wiley & Sons, Inc.) 




(a) 


(b) 


Figure 5.8 

Stress concentration factors for a flat bar with a notch subjected to (a) bending or (b) axial 
load. (From ref. 2; adapted with permission from John Wiley & Sons, Inc.) 



0 0.1 0.2 0.3 0.4 0.5 0.6 

d/b 
(a) 



d/b 

(b) 


Figure 5.9 

Stress concentration factors for a flat plate with a central hole subjected to (a) bending or (b) 
axial load. (From ref. 2; adapted with permission from John Wiley & Sons, Inc.) 
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Figure 5.10 

Theoretical stress concen¬ 
tration factor K t for a ten¬ 
sion-side gear tooth fillet. 
Pressure angle is <p = 20°. 
(From ref. 2; adapted 
with permission from 
John Wiley & Sons, Inc.) 



Figure 5.11 10 

Stress concentration factors K t for a 
crankshaft in bending. Note: When the 
inside of the crankpin and the outside 9 

of the journal are in line, s = 0 When 
the crankpin is closer, s is positive (as 
shown). When the crankpin’s inner 8 

surface is farther away than d/2, S is 
negative. (From ref. 2; adapted with 
permission from John Wiley & Sons, V 

Inc.) 
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p 



p 

& 

h 

d s 

K t 


interfacial pressure 
shaft bending stress 
hub length 
shaft diameter 
theoretical stress 
concentration factor 


Figure 5 .12 

End-of-hub stress concentration fac¬ 
tors for press fit assemblies subjected 
to bending moments. (From ref. 5; 
reprinted by permission of Pearson 
Education, Inc., Upper Saddle River, 
NJ.) 




Figure 5.13 

Stress concentration ef¬ 
fects due to superimposed 
multiple notches. (After 
ref. 6; adapted with per¬ 
mission of The McGraw- 
Hill Companies.) 


the stress concentration factor K tl for the small notch acting alone, it may be multiplied 
times a' n to obtain the actual stress at the small notch root, or 


Next, utilizing (5-27), 


®actual n 


®actual K tl K t2 (Tn 


(5-28) 

(5-29) 


Thus, the combined theoretical stress concentration factor K tc for the multiple notch 
is the product of the stress concentration factors for the two notches considered individu¬ 
ally, giving 

K tc = K tl K t2 (5-30) 

This has been verified photoelastically. 9 The combined fatigue stress concentration factor 
K tc , as for any other stress concentration factor, depends 10 on q and can be calculated from 
(5-92) by substituting K tc for K t in that equation. One technique for estimating a conser¬ 
vative value for K tc is sketched in Figure 5.13(b). The technique assumes the notch of 
5.13(a) to be filled in as shown by the crosshatched area. This leaves a single, deep, nar¬ 
row notch, for which the theoretical stress concentration factor will always be greater than 
the stress concentration factor for the multiple notch. 


Example 5.4 Stress Concentration of Flat Bar 

Two flat bars with different configurations are shown in Figure E5.4 (a) and (b). Each car¬ 
ries a steady load of 4 kN. The material is gray cast iron ASTM A-48 (class 40). Determine 
the maximum stress induced for each of the conditions shown. 


9 Ibid. 10 See (5-91). 
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Example 5.4 
Continues 




30 mm 




AT 5 mmr 


(a) (b) 

Figure ESA 

Flat bars supporting an axial load of 4 kN. 


Solution 

(a) The actual stress at the critical location adjacent to the hole is 

^ (U 

The nominal stress is calculated as 

P 4000 


&nnm : 


A net 0.005(0.040 - 0.010) 


= 26.7 MPa 


From Figure 4.22(b) with d/b = 10/40 = Va 


K t = 2.37 


Hence, the actual stress becomes 

d act = 2.37(26.7) = 63.3 MPa 


(b) The nominal stress is 


P _ 4000 

A net ~ 0.005(0.03) 


26.7 MPa 


From Figure 4.20(b) with r/h = 5/30 = 0.1667, and H/h = 40/30 = 1.333 


K t = 1.7 


Hence, the actual stress becomes 

c r act = 1.7(26.7) = 45.4 MPa 


5.4 Combined Stress Theories of Failure 


Predicting failure, or establishing a combination of material and geometry that will avert 
failure, is a relatively simple matter if the machine part is subjected to a static uniaxial 
state of stress. It is necessary only to have available the simple uniaxial stress-strain curve 
for the material of interest, which can be readily obtained from one or a few simple ten¬ 
sion and compression experiments. For example, if yielding has been established as the 
governing failure mode for a uniaxially stressed machine part, failure of the part would be 
predicted when its maximum normal stress equals or exceeds the uniaxial yield strength 
of the material. 
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If a machine part is subjected to a multiaxial state of stress, the accurate prediction of 
failure becomes far more difficult. No longer can one accurately predict yielding, for ex¬ 
ample, when the maximum normal stress exceeds the uniaxial yield strength because the 
other two principal normal stress components at the critical point may also influence yield¬ 
ing behavior. Further, multiaxial yield strengths are generally unavailable because of the 
time and expense required to determine them experimentally. Therefore, when one desires 
to predict failure, or to pick a combination of material and geometry to avert failure, when 
the machine part is subjected to a multiaxial state of stress, it is usual to utilize an experi- 
mentaly validated theory that relates failure in the multiaxial state of stress to failure by the 
same mode in a simple uniaxial stress test. All such failure prediction theories are based on 
well-chosen loading severity parameters such as stress, strain, or strain energy density. 
Loading severity parameters must be readily calculable in the multiaxial state of stress and 
readily measurable in a simple uniaxial stress test. These theories, called combined stress 
theories of failure, all share a common postulate, namely, that failure is predicted to occur 
when the maximum value of the selected loading severity parameter in the multiaxial state 
of stress becomes equal to or exceeds the value of the same loading severity parameter that 
produces failure in a simple uniaxial stress test using a specimen of the same material. 

Many combined stress failure theories have been proposed, but three have found wide 
acceptance because of their relatively good agreement with experimental results and rea¬ 
sonable simplicity in application. They are 

1. Maximum normal stress theory 

2. Maximum shearing stress theory 

3. Distortion energy theory 


Maximum Normal Stress Theory (Rankine's Theory) 

In words, the maximum normal stress theory may be expressed as: 

Failure is predicted to occur in the multiaxial state of stress when the max¬ 
imum principal normal stress becomes equal to or exceeds the maximum (5-31) 
normal stress at the time of failure in a simple uniaxial stress test using a 
specimen of the same material. 

For a multiaxial state of stress, the maximum principal normal stress is the largest of the 
three roots to the stress cubic equation (5-1), that is, the largest of cr x , ct 2 , and a 3 . The max¬ 
imum normal stress at the time of failure is equal to the uniaxial strength of the material 
corresponding to the governing failure mode. It should be noted that for some materials the 
failure strength under tensile loading may be different from the failure strength under com¬ 
pressive loading. 

With these factors in mind, the word statement (5-31) may be expressed mathemati¬ 
cally as failure is predicted by the maximum normal stress theory to occur if (FIPTOI) 

® 1 “ O'fail-t 

^ Cfail-t (5-32) 

°3 ~ Vfail-t 

or if 

1 — Vfail-c 
— Cfail-c 
~ Vfail-c 


(5-33) 
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Figure 5.14 

Comparison of biaxial ultimate strength data 
for brittle materials with the maximum normal 
stress theory. 


°2 i 


Brittle materials 



<x 3 = 0 
□ Brass 
o Cast iron 
a Cast iron 


°"ult 

- Maximum normal 

stress theory 


where (Jf ai i_ t is the uniaxial tensile ( + ) failure strength of the material and oy ai7 _ c is the uni¬ 
axial compressive (—) failure strength of the material, corresponding to the governing fail¬ 
ure mode (usually yielding or ultimate rupture if loading is static). 

The maximum normal stress theory provides good results for brittle materials, as il¬ 
lustrated in Figure 5.14, but should not be used for ductile materials. 

Maximum Shearing Stress Theory (Tresca-Guest Theory) 

In words, the maximum shearing stress theory may be expressed as 

Failure is predicted to occur in the multiaxial state of stress when the max¬ 
imum shearing stress magnitude becomes equal to or exceeds the maximum (5-34) 
shearing stress magnitude at the time of failure in a simple uniaxial stress 
test using a specimen of the same material. 

For a multiaxial state of stress the maximum shearing stress magnitude is the largest 
of the three principal shearing stresses r ls r 2 , and r 3 , given by (5-2), (5-3), and (5-4). For 
a uniaxial stress test, the only nonzero normal stress component is a principal stress com¬ 
ponent in the direction of the applied force. From (5-2), (5-3), and (5-4), in a uniaxial 
stress test at failure, two of the principal normal stresses are zero and the third is set equal 
to a faih giving 

O'fail 

Tfati = — (5-35) 

With these factors in mind, word statement (5-34) may be expressed mathematically 
as failure is predicted hy the maximum shearing stress theory to occur if (FIPTOI) 

I Til a Wfaii | 

k 2 l 2S Wfaii | (5-36) 

Wi\ - Wfaii I 

where T faU is the largest principal shearing stress at the time of failure in a uniaxial stress 
test, as given by (5-35). 
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Figure 5.15 

Comparison of biaxial yield strength data with 

the maximum shearing stress theory and the 
distortion energy theory . 


The maximum shearing stress theory provides good results for ductile materials, as il¬ 
lustrated in Figure 5.15, but should not be used for brittle materials. 

Distortion Energy Theory (Huber-von Mises-Hencky Theory) 

In words, the distortion energy theory, sometimes called the von Mises theory, may be ex¬ 
pressed as 

Failure is predicted to occur in the multiaxial state of stress when the dis¬ 
tortion energy per unit volume becomes equal to or exceeds the distortion ( 5 . 37 ) 
energy per unit volume at the time of failure in a simple uniaxial stress test 
using a specimen of the same material. 

The distortion energy theory, developed as an improvement over an earlier “total 
strain energy theory,” is based on the postulate that the total strain energy U T stored in a 
volume of stressed material may be divided into two parts: the energy associated solely 
with change in volume, U v , termed dilatation energy, and the energy associated solely 
with change in shape, U d , termed distortion energy. It was further postulated that failure, 
particularly under conditions of ductile behavior, is related only to the distortion energy. 
Thus 


U T = U v + U d (5-38) 

Dividing by the volume in each term, the distortion energy per unit volume may be ex¬ 
pressed as 

u d ~ u T~ u v (5-39) 

An expression for total strain energy per unit volume, u T may be found 11 by calculat¬ 
ing the work done in a triaxial state of stress by forces associated with a u ct 2 , and a 3 , act¬ 
ing over their respective areas, to induce the strains s u s 2 , and e 3 , and their corresponding 
displacements. 

n See equations (4-59), (5-1), (5-15), (5-16), and (5-17). 
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Employing Hooke’s Law, the resulting expression for total strain energy per unit vol¬ 
ume becomes 


U T = ^gL a l + 0-2 + 03 “ 2v(crJ(J 2 + <j 2 o- 3 + o- 3 o-!)] 


(5-40) 


Likewise, the dilatation (volume changing) energy per unit volume may be found 12 as 


3(1 - 2v) 


2 E 


(J \ + (7 2 + (T 3 


(5-41) 


Then substituting (5-40) and (5-41) into (5-39), the distortion energy per unit volume is 
found to be 



1 + v 
3 E 


[(o-i - a 2 ) 2 + (o - 2 - o - 3 ) 2 + (o - 3 - crj) 2 ] 


(5-42) 


To find the distortion energy per unit volume at the time of failure, u d . faih (5-42) is evalu¬ 
ated under uniaxial failure conditions, when two of the principal stresses are equal to zero 
and the third one is equal to the uniaxial failure strength <jf aU . Thus 


Ud-fail 


1 + V 

3 E 



(5-43) 


With (5-42) and (5-43) at hand, word statement (5-37) may be expressed mathemati¬ 
cally as failure is predicted by the distortion energy theory to occur if (FIPTOI) 

t[(o-i - o- 2 ) 2 + (o- 2 - <r 3 ) 2 + (o- 3 - o-i) 2 ] > (r} ail (5-44) 


The distortion energy theory provides very good results for ductile materials, as illus¬ 
trated in Figure 5.15, but generally should not be used for brittle materials. 

The left side of (5-44) is sometimes defined as the square of the von Mises stress , 
the effective stress , or the equivalent uniaxial stress , a eq , giving the expression 


Veq = ~ 1 (02 - 0- 3 ) 2 + ( 0-3 - ff] (5-45) 

For a biaxial state of stress, assuming cr 3 = 0, (5-45) becomes 

a eq = Vaj - o' | a 2 - aj (5-46) 

In terms of rectangular components of stress, the von Mises criterion can be written as 

a eq = - f) 2 + (<Ty - + fz _ O'*) 2 ] + fjxy + fz + f) (5-47) 

and for a biaxial state of stress, assuming that a z = 0 , and t xz = r yz = 0 , we have 


12 For example, see ref. 1, p. 154. 


CT eq = Vo - 2 - <T/T y + O ’ 2 + 


(5-48) 
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Failure Theory Selection 

Evaluation of the three failure theories just discussed in light of experimental evidence 

leads to the following observations: 

1. For isotopic materials that fail by brittle fracture, the maximum normal stress theory 
is the best theory to use. 

2. For isotropic materials that fail by yielding or ductile rupture, the distortion energy 
theory is the best theory to use. 

3. For isotropic materials that fail by yielding or ductile rupture, the maximum shearing 
stress theory is almost as good as the distortion energy theory. 

4. As a rule of thumb , 13 the maximum normal stress theory would be used for isotropic 
brittle materials (materials that exhibit a ductility of less than 5 percent elongation 
in 2 inches) and either the distortion energy theory or maximum shearing stress the¬ 
ory would be used for isotropic ductile materials (materials that exhibit a ductility 
of 5 percent or more in a 2-inch gage length). Where possible, a fracture mechanics 
analysis should be performed. 


Example 5.5 Yielding Failure Prediction Under Static Multiaxial State of Stress 

An aircraft wing flap actuator housing is made of cast magnesium alloy AZ63A-T4 (S u = 
276 MPa, Syp = 97 MPa, e = 12% in 50 mm). The state of stress has been calculated for 
the suspected critical point and is shown in Figure E5.5. Would you predict failure due to 
yielding? 



Figure E5.5 

State of stress at suspected 
critical point. 


Solution 

For the state of stress shown, the stress cubic equation (5-1) becomes 

< t 3 - a 2 (a x + a y + a z ) + <7(a x a y + a y cr z + a z a x - t 2 z ) - (cr x a y a z - (t x t 2 z ) = 0 

Substituting numerical values (in MPa) 

O- 3 - a 2 ( 140 + 70 + 100) + o-[(140)(70) + (70)(100) + (100)(140) - (40) 2 ] 

- [(140)(70)(100) - 140(40) 2 ] = 0 
a 3 - 310a- 2 + 29 200a- - 756000 = 0 


The 5 percent boundary between “brittle” and “ductile” behavior is an arbitrary but widely used rule of thumb. 
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Example 5.5 Since shearing stresses are zero on the x-plane, it is, by definition, a principal plane. Thus, 

Continues 

(j x = (J\ = 140 MPa 

is one of the three solutions of (1). Next, dividing (1) by (<r — 140) gives 

a 2 - 170a + 5400 = 0 


whence 


170 ± V(170) 2 - 4(5400) 

a = -------- = 85 ± 42.7 

2 


giving the solution 


(t 2 = 127.7 MPa and cr 3 = 42.3 MPa 

Since the elongation in 50 mm is given as 12 %, the material may be regarded as ductile; 
and the distortional energy theory of failure, given in (5-7?) is the best theory to use. Thus, 
FIPTOI 

^[(140 - 127.7) 2 + (127.7 - 42.3) 2 + (42.3 - 140) 2 ]>(97) 2 


or FIPTOI 


8495 > 9409 

Since this condition is not met, failure is not predicted at the specified critical point. 


Example 5.6 Thickness of Thin Cylindrical Pressure Vessel 

A thin-walled cylindrical pressure vessel with closed ends is subjected to an internal 
pressure of 100 MPa. The internal diameter is 50 mm. If the yield stress for the pressure 
vessel is 400 MPa, what should be the wall thickness based upon (a) the maximum shear 
stress theory and (b) distortion energy theory? 

Solution 

The stresses in the cylindrical pressure vessel are (see Chapter 9) given as 

_ P4 _pA 

<T ' 2 1 ’ a ‘ At 

where a 1 is the tangential stress and cr x is the longitudinal stress. The tangential and lon¬ 
gitudinal stresses are also the principal stresses in the pressure vessel. Thus, we have 

100(50) 2500 100(50) 1250 

cr i — a t — - — , cr 2 — cr i — ~ — , cr 3 — 0 

2t t 4t t 

(a) Maximum shear stress theory is given as 


Tmax 


cr 3 


2 


2500 - 0 
2 1 


1250 

t 
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Thus, we have 


t = 


1250 

t 

1250 

200 


“yp 


400 


= — =-= 200 


= 6.25 mm 


(b) Distortion energy theory is given by the following: 


<t = Ver? - a x a 2 + oj 


Hence, 


= S yp = 400 H y V2500 2 - 2500 X 1250 + 1250 2 


2165.06 

400 


5.41 mm 


2165.06 

t 


Example 5.7 Maximum Torque in Shaft Subjected to Bending 

A 2.5-in diameter shaft is made of AISI 1020 cold-drawn steel. The shaft is subjected to a 
bending moment of 50,000 in-lb. Determine the maximum torque that can be applied to 
the shaft according to (a) the maximum shear stress theory and (b) the distortion energy 
theory. 

Solution 

The maximum bending stress occurs on the outer surface of the shaft and is given by 


a 


X 


Me 

I 


32 M 
ird 3 


The maximum shear stress due to torque T occurs on the outer surface of the shaft and is 
given by 

_ Tr _ 16 T 

Txy ~J~—p 

7TCI 

The principal stresses are obtained from 


2 

16 


^ 1,2 = -T * \/ 1 V ) + 4 * 


16M 

ird 3 


= —-\m ± Vm 2 + r 2 l 

TTd 31 J 


( 16M\ 2 /16TV 

V 7 rd 3 ) V ird 3 ) 


(a) The maximum shear stress theory states that the maximum shear stress is 


= ^^ = J6 V] ^ VT 2 = ^ 

2 77 d 3 2 
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Example 5.7 
Continues 


Solving for torque T gives 
T = 




- M A 


1/2 


Thus, for AISI 1020 CD S yp = 51,000 psi, we find 


T = 


51,000(77)( 


- 50,000 2 


1/2 


= 60170 in-lb 


(b) The distortion energy theory is given by 


cr = Vcrf - did 2 + 02 = -V4 M 2 + 3 T 


Solving for torque T gives 


7 = 3 


/ 7 rd 3 


7Tu 


4 M 


0 } 


1/2 


For AISI 1020 CD, we find 


-{K 


51,0001 


^( 2 . 5 )- 

16 


- 4(50,000)' 


1/2 


= 69,478 in-lb 


Example 5.8 Factor of Safety for Shaft Subjected to Bending and Torsion 

A simply supported stationary shaft, 50 mm in diameter and made of AISI 1060 hot-rolled 
steel, is subjected to a maximum bending moment of 3000 N-m and a maximum torque of 
2000 N-m. Find the factor of safety corresponding to failure based on the distortion energy 
theory. 

Solution 

The bending stress is given as 


(T x T 3 

1 77 d 

and the maximum torsional stress is 


Me 32 M 32(3000000) 


Txy J 


Tr 16 T 

77 d 3 


7r(50) 3 


16(2000000) 


7t(50) 3 


= 244.46 MPa 


= 81.49 MPa 


The distortion energy theory is given as 

<r eq = Vo- 2 + 3t 2 3, = V(244.46) 2 + 3(81.49) 2 = 282.28 MPa 
For AISI 1060 hot-rolled steel, the yield strength is 372 MPa; thus the factor of safety is 


n d 


372 

282.28 


= 1.32 
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5.5 Brittle Fracture and Crack Propagation; Linear Elastic 
Fracture Mechanics 


When the material behavior of a machine part is brittle rather than ductile, the mechanics 
of the failure process are much different. As described in 2.3, in brittle fracture the part 
separates into two or more pieces due to the breaking of primary interatomic bonds, with 
little or no plastic flow. High-velocity crack propagation from preexisting flaws results in 
sudden and catastrophic failure. If material behavior is clearly brittle, and geometry and 
loading are simple, as for the axially loaded bar of Figure 2.1, failure by brittle fracture 
may be predicted when the stress cr, from (2-1), exceeds the material’s ultimate strength 
S u . Hence failure is predicted to occur if (FIPTOI) 


^ S u (5-49) 

Thus for a uniaxial state of stress the failure prediction expression (5-49) for brittle 
fracture is formally the same as (2-18) for ductile rupture. If the loading is more compli¬ 
cated, and a multiaxial state of stress is produced by the loads, fracture may be predicted 
with reasonable accuracy through the use of the maximum normal stress theory of failure, 
as discussed in 5.4. 

On the other hand, it has now been well established that nominally ductile materials 
may also fail by a brittle fracture response in the presence of cracks or flaws if the combi¬ 
nation of crack size, geometry of the part, temperature, and/or loading rate lies within cer¬ 
tain critical ranges. The prediction of brittle fracture in these circumstances has been based 
on the assumptions that the stress at a crack tip, where failure is initiated, may be calcu¬ 
lated as if the material behavior is linear elastic and the state of stress is two-dimensional; 
thus the procedure is often referred to as linear elastic fracture mechanics (LEFM). 

Three basic types of stress fields have been defined for crack-tip stress analysis, each 
one associated with a distinct mode of crack deformation, as illustrated in Figure 5.16. The 
crack opening mode, Mode I, is associated with local displacement in which the crack sur¬ 
faces move directly apart. Modes II and III are forward sliding and tearing displacements re¬ 
spectively. Mathematical expressions have been developed for the intensity and distribution 
of stress near the crack tips for each of the three modes shown in Figure 5.16. 14 For failure 
prediction purposes, the crack-tip stress-intensity factors, developed from these mathemati¬ 
cal expressions, provide a good measure of the seriousness of loading and geometry in any 
particular case. In general, the expressions for stress intensity factor, K, are of the form 

K = C(j\ // 7tci (5-50) 

where cr is gross-section nominal stress, a is a crack-length parameter, and C is dependent 
upon the type of loading and the geometry away from the crack. Many values of C have 



Figure 5.16 

Basic modes of crack displacement. 

(a) Mode I, (b) Mode II, (c) Mode III. 
(To be read “Mode-one,” “Mode-two,” 
and “Mode-three.”) 


14 See, for example, the Westergaard equations, ref. 1, pp. 54-55. 
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been published, 15 and several typical charts for selecting proper C values for through-the- 
thickness cracks are given in Figures 5.17 through 5.21. Figure 5.22 gives a chart for sur¬ 
face flaw shape parameter Q to be used in finding stress-intensity factors for part-through 
thumbnail-shaped surface cracks. The stress-intensity factor, K , calculated from (5-50), is 
a single-parameter measure of the seriousness of the stress held around the crack tip. The 
magnitude of K associated with the onset of rapid crack extension (initiation of brittle frac¬ 
ture) has been designated as critical stress intensity , K c . Thus failure by brittle fracture may 
be predicted to occur for through-the-thickness cracks if (FIPTOI) 


K = CaV^a > K c 


(5-51) 


or, for thumbnail-shaped surface cracks if 


K = 


1.12 

vi 


d\fim > K c 


(5-52) 


For a given cracked plate, for example the case shown in Figure 5.19, the stress-inten¬ 
sity factor K increases proportionally with gross nominal stress <r, and also is a function of 


Figure 5.17 

Stress-intensity factors K h K n and 
K IU , for center-cracked test speci¬ 
men. (Source: ref. 11, Del Research 
Corp.) 


ttttttt 
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K m - c \u T e 
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4 


— = 0.637 
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alb 


15 See, for example, ref. 7. 










































Brittle Fracture and Crack Propagation; Linear Elastic Fracture Mechanics 235 


ttttttt 



mini 

<j 

K] = CjO" J 7 ,m K u = CjjT ^ttci 



Figure 5.18 

Stress-intensity factors Kj , K n and 
K IU , for double-edge notch test spec¬ 
imen. (Source: ref. 7, Del Research 
Corp.) 
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Figure 5.79 

Stress-intensity factors Kj , for 
single-edge notch test specimen. 
(Source: ref. 7, Del Research 
Corp.) 
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Figure 5.20 

Stress-intensity factors K h 
for single through-the- 
thickness edge crack un¬ 
der pure bending 
moment. (Source: ref. 7, 
Del Research Corp.) 
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Figure 5.21 

Stress-intensity factors K h for a through-the-thickness crack emanating from a circular hole in an infinite plate under 
biaxial tension. (Source: ref. 7, Del Research Corp.) 


Figure 5.22 

Surface flaw shape parameter. (From 
adapted by permission of Pearson 
Education, Inc., Upper Saddle River, 
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Figure 5.23 

Coordinate system for infinite plate containing a 
through-the-thickness crack of length 2a. 


instantaneous crack length a. For the infinite plate shown in Figure 5.23, the crack is ori¬ 
ented so that its plane is perpendicular to the direction of applied uniaxial stress, <r, hence 
pure Mode 1 loading exists. The crack front extends uniformly all the way through the plate, 
so the crack is a through-the-thickness crack. The usual coordinate system for defining state 
of stress at the crack tip is shown. In studying material behavior it has been found that for 
a given material, depending upon the state of stress at the crack tip, the critical stress inten¬ 
sity K c decreases to a lower limiting value as the state of strain at the crack tip approaches 
the conditions of plane strain. 16 This lower limiting value defines a basic material property, 
the plane strain fracture toughness , designated 17 K Ic . Standard test methods have been es¬ 
tablished for the determination of K Ic values. 18 A few data are given in Table 5.2. For the 
plane strain fracture toughness K Ic to be a valid failure prediction criterion for a machine 
part, plane strain conditions 19 must exist at the crack tip; that is, the material must be thick 
enough to ensure plane strain conditions. It has been estimated empirically that for plane 
strain conditions the minimum material thickness B must be 


B > 2.5 



(5-53) 


If the material is not thick enough to meet the criterion of (5-53), plane stress better 
characterizes the state of stress at the crack tip, and K c , the critical stress-intensity factor 
for failure prediction under plane stress conditions, may be estimated using a semiempiri- 
cal relationship for K c as a function of plane strain fracture toughness K Ic and thickness 
B 20 This relationship is 


K r = K h 


L4 | 

(K lc \ 1 

1 H- ~ 

. — 

l B 2 

\ S yp J J 


1/2 


(5-54) 


16 For thick plates the surrounding material constrains the crack-tip zone to near-zero strain in the thickness 
direction, resulting in plane (biaxial) strain. See also 5.2. 

17 To be read “K-one-c.” 

18 See ref. 8. 

19 Plane strain (biaxial state of strain) occurs when nonzero strain components exist in only two coordinate 
directions. See also 5.2. 

20 See refs. 13 and 14. 
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TABLE 5.2 Yield Strength and Plane Strain Fracture Toughness Data for Selected Engineering Alloys 1 


Alloy 

Form 

Test 

Temperature 

S yp 



K Ic 

°F 

°C 

ksi 

MPa 

ksiVin 

MPaVm 

AISI 1045 steel 

Plate 

25 

-4 

39 

269 

46 

50 

AISI 1045 steel 

Plate 

0 

-18 

40 

276 

46 

50 

4340 steel (500°F temper) 

Plate 

70 

21 

217-238 

1495-1640 

45-57 

50-63 

4340 steel (800°F temper) 

Forged 

70 

21 

197-211 

1360-1455 

72-83 

79-91 

D6AC steel (1000°F temper) 

Plate 

70 

21 

217 

1495 

93 

102 

D6AC steel (1000°F temper) 

Plate 

-65 

-54 

228 

1570 

56 

62 

18 Ni maraging steel (300) 

Plate 

600 

316 

236 

1627 

80 

87 

18 Ni maraging steel (300) 

Plate 

70 

21 

280 

1931 

68 

74 

18 Ni maraging steel (300) 

Plate 

-100 

-73 

305 

2103 

42 

46 

A 538 steel 

— 

— 

— 

250 

1722 

100 

111 

2014-T6 aluminum 

Forged 

75 

24 

64 

440 

28 

31 

2024-T351 aluminum 

Plate 

80 

27 

54-56 

370-385 

28-40 

31-44 

6061-T651 aluminum 

Plate 

70 

21 

43 

296 

26 

28 

6061-T651 aluminum 

Plate 

-112 

-80 

45 

310 

30 

33 

7075-T6 aluminum 

— 

— 

— 

75 

517 

26 

28 

7075-T651 aluminum 

Plate 

70 

21 

75-81 

515-560 

25-28 

27-31 

7075-T7351 aluminum 

Plate 

70 

21 

58-66 

400-455 

28-32 

31-35 

Ti-6A1-4V titanium 

Plate 

74 

23 

119 

820 

96 

106 


Trom refs. 10-12. 


As long as the crack-tip plastic zone is in the regime of small-scale yielding, 21 this esti¬ 
mation procedure provides a good design approach. If the plastic zone size ahead of the 
crack tip becomes so large that the small-scale yielding condition is no longer satisfied, an 
appropriate elastic-plastic fracture mechanics (EPFM) procedure would give better re¬ 
sults. For example, a failure assessment diagram 22 might be utilized; however such EPFM 
procedures are beyond the scope of this text. 

To utilize (5-52) as a design or failure prediction tool, the stress-intensity factor K 
must be determined for the particular loading and geometry of the part or structure under 
consideration. The critical stress intensity is set equal to K Ic if the minimum thickness cri¬ 
terion (5-53) is met, otherwise K c is estimated from (5-54). If failure is predicted to occur 
using (5-52), redesign becomes necessary. 


Example 5.9 Brittle Fracture 

Two “identical” support straps of forged 2014-T6 aluminum, shown in Figure E5.8, have 
been inspected and found to contain through-the-thickness cracks. While the total crack 
length of 5 mm is the same for both straps, one strap (case A) involves two edge cracks, 
each with 2.5-mm length, opposite each other; while the other strap (case B) involves a sin¬ 
gle crack of 5-mm length at the center. The straps are of rectangular cross section 50 mm 


2l Small-scale yielding means that the crack-tip plastic zone size is small compared to the dimensions of the 
crack. 

22 See ref. 1, pp. 70-76. 
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660 mm 




\ / 

12.5 mm / 

p 


2.5 mm J 
t 

f 



i 

50 mm 



2.5 mm | 

1 


Case A 


Figure £5.9 

Cracked support bracket. 


/ 


12.5 mm / 



5 mm 

\ 

50 mm 

1 



Case B 


wide and 12.5 mm thick, and the straps are 660 mm in length. The straps are loaded by a 

direct tensile force P = 220 kN in the 660 mm direction, as shown. 

a. For each case, what is the governing failure mode? 

b. For each case, is failure predicted to occur? 

c. Is center cracking or edge cracking the more serious configuration? 

d. How would the analysis change if the strap thickness for each case had been 11 mm 
instead of 12.5 mm? 

Solution 

a. Both yielding and brittle fracture should be checked for each case. Property values 
from Table 5.2 are 

S w = 440 MPa and K lc = 31 MPaVm 

b. For each case, failure by both yielding and brittle fracture must be checked. For case A 
(edge cracks), to check for yielding, from (2-17) FIPTOI 


Where 


P 220,000 

a = -= -= 391 MPa 

A net (0.05 - 0.005)(0.0125) 

Hence, FIPTOI 391 ^ 440. Therefore, failure by yielding is not predicted for case A. 
To check brittle fracture, from (5-52) FIPTOI 

C A a\fim > K c 

Using the Mode I curve from Figure 5.18, with a = 2.5, b = 25, and a/b = 0.1, 
C A V 1 - 0.1 = 1.07 -*• C A = 1.13 


The cross-sectional stress is 


a 



220,000 

(0.05)(0.0125) 


= 352 MPa 
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Example 5.9 Checking for plane strain, from (5-50) plane strain exists if 

Continues 

sa25 (S) j=25 (^) 2=00124 m 

Since B = 12.5 mm, the plain strain condition is met and plain strain conditions pre¬ 
vail. Therefore, 

K c = K Ic = 31 MPaVm 

and the failure criteria may be evaluated as FIPTOI 

CaCtVttu > K c : (1.13)(352)VV(0.0025) s 31 


or if 


35.2 > 31 

Hence, failure by brittle fracture is predicted for case A. 

For case B (center crack), the calculation for yielding failure is identical to the 
case A calculation, since the net cross-sectional area is the same in both cases. 
Therefore, failure by yielding is not predicted for case B. To check brittle frac¬ 
ture, from (5-52) FIPTOI 

C b o-\/ttci > K c 

Using the Mode I curve from Figure 5.17, with a = 2.5, b = 25, and a/b = 0.1, 
C B Vl - 0.1 = 0.96 -» C B = 1.01 

The cross-sectional stress is the same as for case A. Thus, a = 352 MPa. 
Since the material properties are the same as for case A, plain strain conditions 
prevail if and 

K c = K Ic = 31 MPaVm 

Now the failure criteria may be evaluated as FIPTOI 
(1.01)(352)VV(0.0025) > 31 


or if 


31.5 > 31 

Hence, failure by brittle fracture is also predicted for case B, but only by a slight 
margin. 

c. Since the stress-intensity factor K for case A (35.2 MPaV^m) is larger than for case B 
(31.5 MPaV^), the edge-crack case (case A) is more serious. 

d. If the strap thickness is reduced from 12.5 mm to 11 mm, the plain strain condition is 
no longer met and the nominal stress levels increase. To check for yielding for both 
cases A and B we compute the stress 

P 220 000 

(j =-as-= 444 MPa 

A net (0.05 - 0.005)(0.011) 

FIPTOI 444 > 440. Therefore, failure by yielding is predicted for both cases. 
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To check for brittle fracture, since only the thickness has changed, the values 
C A = 1.13 and C B = 1.01 remain unchanged. The gross-section stress for both 
cases becomes 


<r 


220,000 

(0.05)(0.011) 


= 400 MPa 


Since plain strain conditions are not satisfied, from (5-51) 


K c = K Ic 1 




411/2 


31 


1.4 / 31 

(0.011) 2 V 440 


411/2 


35.1 MPa 


Then for case A, FIPTOI 

(1.13)(400) V77(0.0025) > 35.1 


or if 


40.1 > 35.1 

And failure by brittle fracture is again predicted for case A. 
Similarly, for case B, FIPTOI 

(1.01)(400) V77(0.0025) > 35.1 


or if 


35.8 > 35.1 


So brittle fracture is again predicted for case B. 


If loading conditions are more complicated, fracture mechanics procedures are avail¬ 
able but much more involved. For example, in cases where the plane of a through-the- 
thickness crack is oriented at some other angle than perpendicular to the direction of cr, the 
applied stress field induces a combination of Mode I and Mode II loading on the crack. 
Mixed-mode fracture behavior presents an analytical challenge, and a designer would usu¬ 
ally consult with a fracture mechanics specialist to evaluate such a problem. Likewise, a 
specialist would often be involved in addressing problems associated with high-rate load¬ 
ing, cyclic loading, or cases where LEFM conditions are not valid. 


5.6 Fluctuating Loads, Cumulative Damage, and Fatigue Life 


In modern engineering practice, repeated loads, fluctuating loads, and rapidly applied 
loads are far more common than static or quasistatic loads. By far, the majority of engi¬ 
neering design environments involve machine parts subjected to fluctuating or cyclic 
loads. Such loads induce fluctuating or cyclic stresses that often result in failure by fatigue. 
Fatigue is a progressive failure process that involves the initiation and propagation of 
cracks until one reaches an unstable size, triggering a sudden catastrophic separation of 
the affected part into two or more pieces. It is difficult to detect the progressive changes in 
material properties that occur during fatigue stressing, and fatigue failure may therefore 
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occur with little or no warning. Periods of rest, with the fatigue stress removed, do not lead 
to any significant healing or recovery from the effects of the prior cyclic stressing. Hence, 
the damage done during the fatigue process is cumulative. Fatigue failures have been rec¬ 
ognized for over 150 years, but only with the advent of high-speed high-performance ma¬ 
chinery, and the development of the aerospace industry, has widespread attention been 
directed toward trying to understand the fatigue process. 

In recent years it has been recognized that the fatigue failure process involves three 
phases. A crack initiation phase occurs first, followed by a crack propagation phase; finally, 
when the crack reaches a critical size, the terminal phase of unstable rapid crack growth to 
fracture completes the failure process. Traditionally, the models for analysis and prediction 
of fatigue failure have lumped all three phases together in the stress-life (S-N) approach. 
Numerous analytical/empirical procedures and a large database have been developed to sup¬ 
port the S-N approach. More recently, the separate modeling of each phase has been under 
intense development, and prediction models have now been set forth for each phase sepa¬ 
rately. This methodology may be referred to as the fracture-mechanics (F-M) approach. 

Additionally, two domains of cyclic loading have been identified on the basis of 
whether the induced cyclic strains are predominantly elastic or predominantly plastic. 
When cyclic loads are relatively low, strain cycles are confined largely to the elastic range, 
and long lives or high numbers of cycles to failure are exhibited, the domain is called high- 
cycle fatigue. When cyclic loads are relatively high, significant levels of plastic strain are 
induced during each cycle, and short lives or low numbers of cycles to failure are exhibited, 
the domain is called low-cycle fatigue, or cyclic strain-controlled fatigue. Occasionally, a 
machine part may be subjected to intermixed loadings from both domains. High-cycle fa¬ 
tigue predominates in most design environments; therefore high-cycle fatigue analysis will 
be emphasized in this text. Low-cycle fatigue analysis is widely discussed in the literature. 23 


Fluctuating Loads and Stresses 

Fluctuating loads and loading spectra , producing associated stress spectra in a machine 
part, reflect the design configuration and operational use of the machine. Perhaps the sim¬ 
plest fatigue stress spectrum to which a machine element may be subjected is a zero-mean 
sinusoidal stress-time pattern of constant amplitude and fixed frequency, applied for a 
specified number of cycles. Such a stress-time pattern, often referred to as a completely re¬ 
versed or zero-mean cyclic stress, is illustrated in Figure 5.24(a). Using the sketches of 
Figure 5.24, we may define several useful terms and symbols; these include 


®max 
®min 


a 


m 


Va 

Act 

R 

A 


maximum stress in the cycle 
minimum stress in the cycle 

,. , ^ max 4 

mean cyclic stress =- 


alternating stress amplitude = 
range of stress = a max - a min 

J . &min 

stress ratio =- 

& max 

a a 

amplitude ratio = — 


23 For example, see refs. 1 or 15. 
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Figure 5.24 

Several constant-amplitude stress-time pat¬ 
terns of interest, (a) Completely reversed; 

R = — 1. (b) Nonzero-mean stress. 

(c) Released tension; R = 0. 



Any two of the quantities just defined, except the combination a a and Aa or the combi¬ 
nation A and R , are sufficient to completely describe the stress-time pattern. 

A second type of stress-time pattern often encountered is the nonzero-mean spectrum 
shown in Figure 5.24(b). This pattern is very similar to the completely reversed case ex¬ 
cept that the mean stress is either tensile or compressive, in any event different from zero. 
The nonzero-mean case may be thought of as a static stress equal in magnitude to the mean 
a m with a superposed completely reversed cyclic stress of amplitude a a . 

A special case of nonzero-mean stress, illustrated in Figure 5.24(c), is often encoun¬ 
tered in practice. In this special case the minimum stress ranges from zero up to some ten¬ 
sile maximum and then back to zero. This type of stressing is often called released tension. 
For released tension it may be noted that a m = cr max / 2. A similar but less frequently en¬ 
countered stress-time pattern is called released compression, where cr max = 0 and 
erm er m in/ 2 - 

More complicated stress-time patterns are illustrated in Figure 5.25. In Figure 5.25(a) 
the mean stress is zero but there are two (or more) different stress amplitudes mixed to¬ 
gether. In Figure 5.25(b), not only does the stress amplitude vary, but also the magnitude of 



(a) 



(b) 


Figure 5.25 
Stress-time patterns in 
which the amplitude 
changes or both mean 
and amplitude change to 
produce a more compli¬ 
cated stress spectrum. 

(a) Zero mean, changing 
amplitude, (b) Changing 
mean and amplitude. 
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Figure 5.26 

A quasi-random stress-time pattern that 
might be typical of an operational aircraft 
during any given mission. 



the mean stress periodically changes, approaching a more realistic condition. Figure 5.26 il¬ 
lustrates a realistic stress-time pattern, as might be observed, for example, in an airframe 
structural member during a typical mission including refueling, taxi, takeoff, gusts, maneu¬ 
vers, and landing. To effectively predict and prevent fatigue failures, procedures for assess¬ 
ing damage due to all of these various stress-time spectra are required. 

Fatigue Strength and Fatigue Limit 

Designing machine parts or structures that are subjected to fatigue loading is usually based 
on the results of laboratory fatigue tests using small polished specimens of the material of 
interest. To be successful in using small-specimen data, a designer must be aware of the 
merits and limitations of such data, including influences of size, surface finish, geometry, 
environment, speed, and many other factors. Basic fatigue data in the high-cycle life range 
are usually displayed on a plot of cyclic stress versus life. These plots, called S-N curves , 
constitute design information of fundamental importance for machine parts subjected to 
cyclic or repeated loadng. 

Figure 5.27 illustrates the characteristic appearance of an S-N curve for a ferrous ma¬ 
terial. The curve is constructed by fitting it to data collected from laboratory tests of a large 
number of specimens at various cyclic stress amplitudes and zero-mean stress. Because of 


Figure 2.27 too 

Plot of stress-life (S-N) data as they might 
be collected by laboratory fatigue testing 
of a new alloy. 80 
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Figure 5.28 

Distribution of fatigue specimen 
failures at a constant stress level 
as a function of logarithm of 
life. 


the scatter of fatigue life data at any given stress level, the construction of the most ap¬ 
propriate curve through the data becomes a design issue of importance. A statistical de¬ 
scription of fatigue failure data is usually used to facilitate construction of the most 
appropriate S-N curve. Using this approach, one would construct for each test stress level 
a histogram , such as the one shown in Figure 5.28, which shows the distribution of fail¬ 
ures as a function of the logarithms of lives for the sample tested. Computation of the sam¬ 
ple mean and standard deviation 24 permits estimation of the probability of failure (P) at 
each stress level. Points of equal probability of failure may then be connected to obtain 
curves of constant probability of failure on the S-N plot. A family of such S-N-P curves is 
shown in Figure 5.29 for 7075-T6 aluminum alloy. It is also of interest that the reliability 
(R) is defined to be 1 minus the probability of failure; hence R = (1 —P). Thus in Figure 
5.29, the 10 percent probability of failure curve (P = 0.10) may alternatively be designated 
as the 90 percent reliability curve (R = 0.90), and some references may be found in the lit¬ 
erature to these so-called R-S-N curves. 

Usually, references in the literature to “the” S-N curve refer to the mean or 50 per¬ 
cent probability of failure curve unless otherwise specified. The mean S-N curves of 
Figure 5.30 distinguish two commonly observed types of response to cyclic loading. The 



10 4 10 5 10 6 10 7 10 8 
Cycles to failure, N 


Figure 5.29 

Family of S-N-P curves, or R-S-N curves for 7075-T6 
aluminum alloy. Note: P = probability of failure; R = relia¬ 
bility = 1 — P. (Adapted from ref. 16, p. 117, Copyright © 
1969, by permission from John Wiley & Sons, Inc.) 



Figure 5.30 

Two types of material response to cyclic loading. 


24 See, for example, ref. 1, Ch. 9. 
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Figure 5.31 

Effect of material composition on the 
S-N curve. Note that ferrous and titanium 
alloys exhibit well-defined fatigue limits 
whereas other alloy compositions do not. 
(Data from refs. 17 and 18.) 
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Cycles to failure 


ferrous alloys and titanium exhibit a steep branch in the shorter life range, leveling off 
to approach a stress asymptote at longer lives. This stress asymptote is called the fatigue 
limit , Sf or S 'f (formerly called endurance limit) and is the stress level below which an 
infinite number of cycles could theoretically be sustained without failure. The symbol Sf 
is commonly used to denote fatigue limit of an actual machine part , while S 'f is fatigue 
limit of the material based on small polished specimens. The nonferrous alloys do not 
exhibit an asymptote, and the curve of stress versus life continues to drop off indefi¬ 
nitely. For such alloys there is no fatigue limit , and failure as a result of cyclic stress is 
only a matter of applying enough cycles. To characterize the failure response of most 
nonferrous materials, and of ferrous or titanium alloys in the finite life range, the term 
fatigue strength at a specified life , Sjy or S' N , is used; S N denotes finite-life fatigue 
strength of a component , while S' N is finite-life fatigue strength of the material based on 
small polished specimens. The specification of fatigue strength without specifying the 
corresponding life is meaningless. The specification of a fatigue limit always implies in¬ 
finite life. Use of these terms is illustrated in Figure 5.31, which depicts S-N curves for 
several alloys. 


Estimating S-N Curves 

When designing a machine part to be subjected to fluctuating loads it is essential to obtain 
an S-N curve for the candidate material for which heat treatment, operating temperature, 
and other operating conditions match the application. If such an S-N curve can be found in 
the literature or in an available database, the curve may be used directly to design the part. 
If not, it may be possible to find an S-N curve for small polished specimens of the candi¬ 
date material and modify it by using a series of appropriate factors as discussed later in 
connection with Table 5.3. If S-N data cannot be found for the candidate material it be¬ 
comes necessary either to perform laboratory fatigue tests on small polished specimens 
(a lengthy and expensive process) or to attempt to estimate the basic S-N curve for the 
material. Although there is no direct connection between static strength and fatigue 
strength, purely empirical relationships have been documented by examining large bodies 
of existing data. These empirical observations allow a designer to estimate S-N curves 
from static ultimate tensile strength ( S ut ) values. These estimated S-N curves are often 
accurate enough for the design task at hand, especially at the preliminary design stage. For 
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TABLE 5.3 Strength-Influencing Factors That May Affect S-N Curves 


Influencing Factor 

Symbol 

Approximate 

Range 

“Typical” Value 1 

Material composition 

— 

— 

Specific S-N data required 

Heat treatment 

— 

— 

Specific S-N data required 

Operating temperature 

— 

— 

Specific S-N data required 

Grain size and direction 

If 

gr 

0.4-1.0 

1.0 

Welding 

If 

0.3-0.9 

0.8 

Geometrical discontinuity 

k f 

0.2-1.0 

Reciprocal of Kf, see 4.8 

Surface condition 

If 

IV sr 

0.2-0.9 

0.7 

Size effect 

ksz 

0.5-1.0 

0.9 

Residual surface stress 

If 

^rs 

0.5-2.5 

Specific data required; see 4.11 

Fretting 

k fr 

0.1-0.9 

0.35 if fretting exists, 1.0 if no fretting; 
also see 2.11 

Corrosion 

If 

i\ cr 

0.1-1.0 

Specific data required 

Operating speed 

K P 

0.9-1.2 

1.0 

Strength reliability required 

k r 

0.7-1.0 

0.9; also see Table 2.3 

Configuration of stress-time pattern 

— 

— 

See later section titled Cumulative 

Damage Concepts and Cycle Counting 

Nonzero-mean stress 

— 

— 

See later section titled Nonzero-Mean Stress 

Damage accumulation 



See later section titled Cumulative 

Damage Concepts and Cycle Counting 


'These “typical” values may be used for solving problems in this text or for making preliminary design estimates when actual conditions are 
poorly known. However, any critical design situation would require a literature/database search or supporting laboratory experiments to establish 
more accurate values. 


wrought-ferrous alloys, the mean S-N curve may be estimated for polished specimens by 
plotting on a semilog plot of strength versus log life the following: 

1. Plot S' N = S ut at N = 1 cycle. 

2. Plot Sf = 0.5 S ut at N = 10 6 cycles if S ut < 200 ksi 
or, 

Plot Sf = 100 ksi at N = 10 6 cycles if S ut > 200 ksi 

3. Connect (1) with (2) by a straight line. 

4. Construct a horizontal straight line from (2) toward very long lives. 

This estimated mean (50 percent reliability) S-N curve may be supplemented if de¬ 
sired by using the standard normal variable X of Table 5.4 to calculate and construct esti¬ 
mated S-N curves with higher reliabilities. 


TABLE 5.4 Strength Reliability Factors as a Function of Reliability 
Level 


Reliability R 
(percent) 

Corresponding Standard 
Normal Variable X 
(see Table 2.9) 

Strength Reliability 
Factor k r 

90 

1.282 

0.90 

95 

1.645 

0.87 

99 

2.326 

0.81 

99.9 

3.090 

0.75 

99.995 

3.891 

0.69 
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For cast irons and cast steels a similar procedure may be used except that 
S' f = 0.4 S ut at N = 10 6 cycles if S ut < 88 ksi 


or 


S'f = 40 ksi if S ut > 88 ksi 

Guidelines for other alloys have been published as follows: 25 


Titanium alloys: Sf 
Aluminum alloys: S' N 
Magnesium alloys: S' N 

= 0.45S„ ( 

= 0.4 S ut 
= 0.35 S ut 

to 

0.65^ 

at 

at 

at 

N = 10 6 cycles 

N = 5 X 10 8 cycles 
N = 10 8 cycles 

Copper alloys: S' N 

= 0.2 S ut 

to 

0.5S„, 

at 

N = 10 8 cycles 

Nickel alloys: S' N 

= 0.3 S ut 

to 

0.5S h( 

at 

N = 10 8 cycles 


It is very important to emphasize that these guidelines provide only estimates , to be used 
when pertinent S-N data are unavailable, and must be used with caution. 

Stress-Life ( S-N) Approach to Fatigue 

For uniaxial states of cyclic stress the traditional S-N approach to fatigue design is straight¬ 
forward in concept. Data are required for the material of interest, including effects of size, 
surface finish, geometry, and environment, and information is required on the cyclic load¬ 
ing to be applied to the part. From these, the size and shape of the part to provide the de¬ 
sired cyclic life may be determined. To implement the S-N approach, material data are best 
presented as S-N curves, such as those shown in Figure 5.30 or 5.31, or Figure E5.10 of 
Example 5.10. For example, if a part to be made of 2024-T4 aluminum has a design life 
requirement of 5 X 10 8 completely reversed stress cycles, its fatigue strength S' N at a life 
N = 5 X 10 8 cycles may be read from Figure 5.31 as 20,000 psi. This assumes that the 
part operates under conditions matching the conditions used to obtain the data for Figure 
5.31. Also, since S N=5xm * = 20,000 psi is a fatigue failure strength, an appropriate safety 
factor would typically be imposed for safe operation, as discussed in Chapter 2. 

Factors That May Affect S-N Curves 

As already mentioned, published S-N curves, unless otherwise labeled, are typically mean 
value curves for small polished specimens. Actual machine parts subjected to cyclic or fluc¬ 
tuating stress levels exhibit S-N responses different from small, polished-specimen curves, 
depending upon differences in composition, processing, environment, and operational fac¬ 
tors. Thus the fatigue strength or fatigue limit of an actual machine part (S N or Sf) is nearly 
always different (usually lower) from the fatigue strength or fatigue limit S' N or S ^read from 
a published S-N curve for small polished specimens of the same material. Furthermore, pub¬ 
lished S-N curves for small polished specimens are somewhat dependent upon the test 
method used to obtain the underlying data. The test method of choice in modern practice is 
to use computer-controlled closed-loop axial push-pull fatigue testing machines that apply 
uniformly distributed cyclic direct stresses to small polished specimens. Most S-N data pre¬ 
sented in this textbook may be assumed to have been produced by uniform uniaxial cyclic 
stresses on small polished specimens, unless otherwise noted. Older small specimen S-N 
data found in the literature may have been obtained from cyclic bending tests, using either 
rotating or reciprocating bending machines. If such bending data are used as the basis for 

25 See ref. 6 
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obtaining fatigue strength or fatigue limit estimates for an actual machine part, as in (5-55) 
or (5-56), the fatigue limit values from the cycle bending data should be reduced by about 
10 percent (because of the bending-stress gradients; see Figure 4.3) to be compatible with 
the uniformly stressed push-pull small specimen data used in this textbook. 

The factors that may cause differences between the fatigue response of actual machine 
parts and small polished specimens are shown in Table 5.3 and briefly discussed in the fol¬ 
lowing paragraphs. More detailed discussions may be found in the fatigue literature. 26 

To obtain an estimate for fatigue strength or fatigue limit of an actual machine part 
from small-polished-specimen S-N data and a knowledge of operational and environmen¬ 
tal requirements for the part, the usual procedure is to utilize an expression of the form 

S f = kooS' f (5-55) 

or 

Sj\f — k N S' N (5-56) 

depending upon whether fatigue limit, or fatigue strength corresponding to a design life of 
N cycles, is required. From Table 5.3 the expressions for koo and k N may be written as the 
products of pertinent influencing factors. Thus 

koo (kgrk we^f^srksz^rs^frkcrkspky-) oo (5-57) 

and 

^ n ~ ( k gr k we kfk sr k sz k rs kf r k cr k S p K)n (5-58) 

It should be recognized that the various factors may have different values for finite-life de¬ 
sign requirements as compared to infinite-life design specifications. Also, it should be 
noted that these procedures assume that no interactions occur between or among the vari¬ 
ous strength-influencing factors, an assumption not always justified. For example, the sur¬ 
face condition effect and corrosion effect may interact, or other effects may interact to 
produce unknown influences on the fatigue strength or fatigue limit. For this reason it is 
always advisable in practice to subject the final design to full-scale prototype testing un¬ 
der simulated service conditions. 

Approximate ranges and “typical” values for many of the strength-influencing factors 
are shown in Table 5.3. While these may be helpful guidelines in approximating the fa¬ 
tigue strength or fatigue limit of a machine part, more specific values are usually required 
for any critical design situation, and a search of fatigue databases or fatigue literature is 
common practice. The following brief remarks relative to the entries of Table 5.3 are of¬ 
fered as additional guidance. 

Material composition is the most basic factor in determining fatigue strength. As 
noted in the discussion of Figure 5.30, materials divide themselves into two broad groups 
with respect to S-N failure response. The ferrous alloys and titanium exhibit a rather well- 
defined fatigue limit, which is well established by the time 10 7 cycles of stress have been 
applied. The other nonferrous alloys do not exhibit a fatigue limit at all, and their S-N 
curves continue to fall off at lives of 10 8 , 10 9 , and larger numbers of cycles. 

Heat treatment is also a strong influencing factor on the fatigue strength, just as it is 
in the case of static strength. Small-polished-specimen S-N data for which the specimen 
heat treatment corresponds to that proposed for the actual part under consideration is very 
useful in determining fatigue response of the part. 

26 See for example, ref. 1, 6, 17, 19-23. 
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Operating temperature may have a significant influence on fatigue strength. 
Generally, the fatigue strength is somewhat enhanced at temperatures below room tem¬ 
perature and diminished at temperatures above room temperature. In the range of temper¬ 
ature from around zero on the Fahrenheit scale to about one-half the absolute melting 
temperature, the effects of temperature are slight in most cases. At higher temperatures the 
fatigue strength diminishes significantly. Further, alloys that exhibit a fatigue limit at room 
temperature tend to lose this characteristic at elevated temperatures, making infinite-life 
design at elevated temperatures impossible. 

Grain size and grain direction may play a significant role in fatigue strength. Fine¬ 
grained materials tend to exhibit fatigue properties that are superior to coarse-grained mate¬ 
rials of the same composition. This superiority becomes less significant at elevated 
temperatures when the characteristic room temperature transgranular cracking gives way to 
intergranular cracking paths. When the cyclic loading direction is across the grain (transverse 
direction) of an anisotropic specimen or machine part, fatigue strength properties are inferior 
to fatigue strength when the loading direction is along the grain (longitudinal direction). 

Welding produces a metallurgically nonhomogeneous region ranging from unheated 
parent metal through the heat-affected zone (HAZ), to the weld metal zone. In some cases 
the entire welded joint may be post-heat-treated, in which case the structure of weld metal, 
HAZ, and parent metal may become nearly identical; in some critical cases postweld ma¬ 
chining operations may be used to restore geometrical uniformity. Even with such care 
(and expense), welded joints (as well as bolted, riveted, or bonded joints) tend to have a 
fatigue strength inferior to that of a monolithic part of the same material. Factors con¬ 
tributing to fatigue strength reduction in welded joints, in addition to the gradient in ho¬ 
mogeneity across the weld zone, include cracking in weld metal or base metal due to 
postcooling shrinkage stresses, incomplete penetration, lack of fusion between weld metal 
and parent metal from prior welding passes, undercut at the edge of weld metal deposit, 
overlap of weld metal flowing beyond the fusion zone, slag inclusions, porosity, misshapen 
welds, or welds with surface defects. Some of the effects of welding on fatigue strength 
properties are illustrated in Figure 5.32. 


Figure 5.32 

Effects of welding detail on the S-N 
curve of structural steels, with yield 
strengths in the range 30,000-52,000 
psi. Tests were released tension 
(iT m in = 0) (Data from ref. 23.) 
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Stress concentration effects due to geometrical discontinuities , such as changes in 
shape or joint connections, may seriously diminish the fatigue strength of a machine part, 
even if the part is made of a ductile material. The severity of notches, holes, fillets, joints, 
and other stress raisers depends upon the relative dimensions, type of loading, and notch 
sensitivity of the material. A detailed discussion of stress concentration is presented in 5.3. 

Surface condition is a very important factor in fatigue strength of a machine part 
since a very high proportion of all fatigue failures nucleate at the surface. Irregular sur¬ 
faces and rough surfaces generally exhibit inferior fatigue properties as compared to 
smooth or polished surfaces, as illustrated in Figure 5.33. Cladding, plating, or coating 
may reduce the fatigue strength of the plated member because fatigue nuclei initiated in 
the plating continue to propagate into the base metal. This may cause a reduction in fa¬ 
tigue strength, as illustrated in Figure 5.34. Usually, however, the corrosion protection af¬ 
forded by the plating more than offsets the strength loss due to plating. 


Hardness, BHN 
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Figure 5.33 

Reduction of fatigue strength due to 
surface finish (steel parts). (From ref. 6, 
reproduced with permission of The 
McGraw-Hill Companies.) 



Figure 5.34 

Effects of several electrodeposited coatings on the 
S-N curve of low-alloy steel at room temperature 
under axial tension-tension loading with R = 0.02. 
Static ultimate strengths: unplated, 172,300 psi; 
corronized, 177,700 psi; hard nickel plated, 
176,100 psi; electroless nickel plated, 182,100 psi; 
hard aluminum plated, 162,400 psi. (Data from 
ref. 16, Copyright © 1969, adapted by permission 
from John Wiley & Sons, Inc.) 
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Figure 5.35 

Effects of shot-peening and/or presetting on the 
S-N curve of hot wound helical coil springs made of 
0.9 percent carbon steel. Spring details: hardness = 
Vickers DPH 550; wire diameter = 1 / 1 inch; mean coil 
diameter = 2% inches; number of turns = 6; free 
length = 5% 6 or 6 inches. (Data from ref. 24.) 
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Figure 5.36 

Effects of cold-rolling threads before and after heat treat¬ 
ment on the S-N curve for 220,000 psi ultimate strength 
bolts, (a) Rolled after heat treatment, (b) Rolled before heat 
treatment. (Data from ref. 16, Copyright © 1969, by 
permission from John Wiley & Sons, Inc.) 


Larger specimens and machine parts exhibit a size effect , having lower fatigue 
strengths than smaller specimens of the same material. For example, the fatigue strength 
of a 6-inch-diameter machine part might be as much 15 or 20 percent lower than a one- 
half-inch-diameter specimen of the same material. 

Residual stresses in the surface layer, whether induced intentionally or accidentally, may 
play an extremely important role in the overall fatigue response of a specimen or machine 
part. If the induced residual surface stresses are tensile, the fatigue strength is diminished. If 
the residual surface stresses are compressive, the fatigue strength is improved. Three com¬ 
mon methods of inducing compressive residual surface stresses are shot-peening, cold¬ 
rolling, and presetting. 27 Figures 5.35 and 5.36 illustrate these effects. It is also notable that 
surface treatments such as shot-peening and cold-rolling not only improve the mean fatigue 
strength, but reduce the scatter (standard deviation) as well. 

Fretting at the contacting surfaces of joints or connections may lead to a very significant 
reduction in fatigue strength of a machine part. As discussed in 2.9, the prediction of the fret¬ 
ting effect is difficult because many factors are involved. Experimental testing of the assem¬ 
bly under actual service conditions and loading should be undertaken for critical parts. 

Corrosion tends to lower the fatigue strength, often by a large amount. Corrosive ef¬ 
fects are specific to the combination of material composition and operating environment. 
Experimental testing of the machine part under service environment and loading should be 
employed for critical parts. 

27 Shot-peening involves uniform bombardment of the surface with a high-velocity stream of small steel spheres. 
Cold-rolling is accomplished by pressing a hard contoured-roller against the surface to be treated and uniformly 
rolling and translating the roller to cover the area to be cold worked. Presetting involves the application of a 
static overload on the part, in the direction of operational loads, to induce local yielding at stress concentrations 
and subsequent compressive residual stresses at these sites upon release of the static overload. Also, see 4.9. 
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Operating speeds in the range from about 200 cycles per minute (cpm) to about 7,000 
cpm appear to have little effect on fatigue strength. Below 200 cpm there is often a small 
decrease in fatigue strength, and in the range 7,000 cpm to around 60,000 cpm many mate¬ 
rials exhibit a significant increase in fatigue strength. Rest periods have no effect on fatigue 
strength. 

Strength reliability , or reliability of fatigue strength, is usually based on the stresswise 
distribution of S-N-P data for small polished specimens of a particular material. Assuming 
the stresswise distribution of fatigue strengths to be normal is a reasonable and convenient 
procedure for estimating the strength reliability factor. 28 From any table for the cumulative 
distribution function 29 for the standard normal variable X , where 


X = 


S'fW ~ S f 

(T 


(5-59) 


and using the empirical estimate 30 for standard deviation of fatigue strength to be 
<t = 0.08 Sf, the strength reliability factor may be expressed as 

k r = 1 - 0.08X (5-60) 

Based on this relationship, values for k r are shown in Table 5.4 for a range of strength re¬ 
liability values, R. If statistical data were available for the distributions of other factors 
listed in Table 5.3, a similar reliability approach could be taken for these influencing fac¬ 
tors, but such data are not generally available. 

In practice, the configuration of the stress-time pattern may take many forms, includ¬ 
ing sinusoidal, reversed ramp, superposed ripples, or distorted peaks as shown in Figure 
5.37. Generally, the fatigue failure response seems to be relatively insensitive to changes 






Figure 5.37 
A variety of stress-time 
patterns used in evaluat¬ 
ing fatigue behavior. 

(a) Completely reversed 
sinusoid, (b) Completely 
reversed ramp. 

(c) Superposed ripples. 

(d) Secondary peaks. 

(e) Distorted peaks. 


28 Although the Weibull distribution is widely agreed to be more accurate, its use is more cumbersome. 
See ref. 1, Ch. 9, for example. 

29 See, for example, Chapter 2, Table 2.9. 30 See ref. 24. 
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Figure E5.10 

Estimated S-N curves for 
fatigue strength reliability 
levels of 50 percent (mean) 
and 99.9 percent. 


in wave shape as long as peak value and period are the same. When practical, however, it 
is better to account for each and every stress reversal by using a good cycle counting 
method such as the rain flow method discussed in the section “Cumulative Damage 
Concepts and Cycle Counting” below. 

The magnitude of any nonzero-mean stress has an important influence on fatigue 
response and is discussed in detail in the next section. Likewise, the accumulation of 
fatigue damage caused by cyclic loading, or cumulative damage , is fully discussed in a 
later section. 


Example 5.10 Estimating Fatigue Properties 

A wrought carbon-steel alloy is known to have the static properties S u = 76,000 psi, S yp = 
42,000 psi, and e (2 inches) =18 percent, but fatigue properties cannot be located for the 
material. It is necessary to quickly estimate the fatigue properties for the preliminary design 
of a machine part for which the fluctuating loads will induce a stress spectrum with cyclic 
amplitudes in both the finite-life range and the infinite-life range. 

a. How could the basic small-polished-specimen “ mean S-N curve” be estimated for this 
material? 

b. How could the R = 99.9 percent reliability S-N curve be estimated for the material? 

Solution 

a. Using the estimation guidelines for wrought-ferrous alloys with Su, < 200 ksi (since 
S ut = S u = 76,000 psi), 

S}= 0.5 (S u ) = 0.5(76,000) = 38,000 psi 
5J^ = i = S u = 76,000 psi 

Plotting these values on semilog coordinates and connecting them in accordance 
with the guideline procedures gives the mean (.R = 50 percent) S-N curve shown in 
Figure E5.10. 

b. To obtain an estimate of the R = 99.9 percent reliability S-N curve, the expression for 
standard normal variable X given in (5-59) may be utilized. Since the R = 99.9 percent 
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curve must lie below the mean curve, 5} (. R = 99.9) is less than the mean value Sf. 
Hence, for this case, from (5-59) 

S' f - S'AR = 99.9) 

X* - x - ( 1 ) 

(7 

Therefore 


Sf = 38,000 psi 

and using the estimated standard deviation a discussed just after (5-59), 
a = 0.085}= 0.08(38,000) = 3040 psi 
with the Table 5.4 value for X corresponding to R — 99.9 percent reliability 

X(R = 99.9) = 3.09 

equation (1) may be solved for 5} (R = 99.9) to give 

Sf(R = 99.9) m S' f - Xa 
or 

S' f (R = 99.9) = 38,000 - 3.09(3040) = 28,600 psi 

Plotting Sf(R = 99.9) = 25,600 psi on the S-N coordinates of Figure E5.10 at 10 6 cycles, 
and connecting to the point $ut at N = 1 cycle, results in the estimated 99.9 percent relia¬ 
bility S-N curve, as shown. 


Example 5.11 Failure Mode Assessment 

An axially loaded straight cylindrical bar of diameter d = 12.5 mm is to be made of 2024- 
T4 aluminum with ultimate strength of S u = 469 MPa, yield strength S yp = 331 MPa, and 
fatigue properties shown in Figure 5.31. The bar is to be subjected to a completely reversed 
axial force of 27 kN, and must last for at least 10 7 cycles. 

a. What is the governing failure mode? 

b. Is failure predicted to occur? 

Solution 

a. The two most probable candidates for governing failure mode are yielding and fatigue. 
Both should be calculated to determine which one governs. 

b. For yielding, from (2-17) FIPTOI 


Now 


F} _ 4 F 
A 0 77 d 2 


4(27000) 

7r(0.0125) 2 


= 220 MPa 
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Example 5.11 
Continues 


Hence, for yielding, FIPTOI 

220 > 331 

Therefore, failure by yielding is not predicted to occur. 

For fatigue, at a design life requirement of N d = 10 7 cycles, FIPTOI 

(T max — Sat=io 7 

From Figure 5.31, using the curve for 2024-T4 aluminum for 10 7 cycles, the fatigue 
strength at N = 10 7 cycles may be read as 


S N=1 o? = 23 ksi » 159 MPa 

Hence, for fatigue at a design life requirement of N d = 10 7 cycles, FIPTOI 

220 > 159 

Therefore, failure by fatigue is predicted. That is, fatigue failure would be expected before 
the design life of 10 7 cycles is achieved. 


Example 5.12 Design for Infinite Fatigue Life 

An axially loaded straight cylindrical bar is to be made of 1020 steel, with fatigue proper¬ 
ties as shown in Figure 5.31. The bar is to be subjected to a completely reversed axial force 
of 7000 lb maximum. Fatigue is the governing failure mode. 

If infinite life is desired for this part, what is the minimum diameter that the bar should 
be made? 

Solution 

Since fatigue is the governing failure mode and infinite life is desired, the stress in the bar 
must be just below the fatigue limit, read from Figure 5.31 for 1020 steel as 

Sf = 33,000 psi 

The maximum stress in the bar is 


^ max 




or 


d = 


4 f 

' 1 m 


Setting cr max = Sf = 33,000 psi (incipient failure), 


/ 4(7000) 

d = ,/ & 0.52 inch 

V 77(33,000) 

Thus the minimum diameter for infinite life would be 0.52 inch. Again, it must be 
cautioned that this assumes that the S-N curve of Figure 5.31 properly reflects the con¬ 
ditions of the axially loaded bar in terms of size, surface finish, geometry, environment, 
and other influence factors. Further, statistical scatter or reliability level should also be 
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considered in establishing the bar diameter, using techniques such as the ones described 
in Example 5.9(b). As a practical matter, the small-polished-specimen S-N curve in 
Figure 5.31 would typically be modified by using appropriate factors from Table 5.3 to 
reflect the operating conditions, and a safety factor would be used to protect against re¬ 
maining uncertainties. 31 


I Example 5.13 Estimating Fatigue Properties of a Part 

A wrought-carbon steel alloy with S u = 524 MPa, S yp = 290 MPa, and <?(50 mm) = 18%, 
is to be considered a candidate material for a proposed machine part to be used under the 
following conditions: 

a. The part is to be lathe-turned from a bar of the wrought-steel alloy. 

b. The part is uniform in shape at the critical point. 

c. Operating speed is 3600 rpm. 

d. A very long life is desired. 

e. A strength reliability factor of 99.9% is desired. 

It is desired to make a preliminary design estimate of the pertinent fatigue properties of the 
candidate material so that approximate size of the part can be estimated. 

Solution 

Since a very long life is desired, the fatigue property of primary interest is the fatigue limit 
of the part, Sf. 

From (5-52) and (5-54), 

S f = kooSf and k oo = (k gr k we kjk ST k sz k rs k fr k cr k sp k r ) 

Using estimation guidelines for wrought-ferrous alloys with S u < 200 ksi(1379 MPa) 
S f = 0.5(SJ = 0.5(524) = 262 MPa 

To evaluate koo, each of the influence factors must be evaluated. Referring to Table 5.3 
and the discussion following the table: 

k gr =1.0 (from Table 5.3) 

k we =1.0 (no welding anticipated) 

kf= 1.0 (uniform shape specified) 

k sr = 0.65 (see Figure 5.33) 

k sz = 0.9 (size unknown; use Table 5.3) 

k rs =1.0 (no information available; later review essential) 

kf r = 1.0 (no fretting anticipated) 

k cr =1.0 (no information available; later review essential) 
k sp =1.0 (moderate; use Table 5.3) 
k r = 0.75(from Table 5.4 for R = 99.9) 

Now we evaluate ^ as 

k^ = (1.0) (1.0) (1.0) (0.65) (0.9) (1.0) (1.0) (1.0) (1.0) (0.75) = 0.43875 = 0.44 
31 See Chapter 2. 
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Example 5.13 The fatigue limit is therefore 

Continues 

S f = 0.44(262) = 115.3 MPa 

Thus, the estimated fatigue limit for the part is 115.3 MPa. Using an appropriate 
safety factor with Sf would allow calculation of the approximate size of the part to produce 
infinite life. 


Nonzero-Mean Stress 

Fatigue failure data collected in the laboratory for small polished specimens are most of¬ 
ten obtained for completely reversed, or zero-mean alternating stresses. Many, if not most, 
service applications involve nonzero-mean cyclic stresses. It is very important, therefore, 
to be able to predict the influence of mean stress on fatigue behavior so that completely re¬ 
versed small-polished-specimen S-N data may be utilized to design machine parts sub¬ 
jected to nonzero-mean cyclic stresses. 

High-cycle fatigue data collected from a series of experiments devised to investigate 
combinations of alternating stress amplitude cr a and mean stress a m may be characterized 
by a plot of a a versus a m for any specified failure life of N cycles, such as Figure 5.38. As 
shown, the failure data points typically tend to cluster about a curve that passes through the 
point cr a = S N at cr m =0 and the point a m = S u at a a = 0. As shown, failure is very sensi¬ 
tive to the magnitude of mean stress in the tensile mean-stress region but rather insensitive 
to the magnitude of mean stress in the compressive mean-stress region. Nonzero-mean fail¬ 
ure data are available in the literature for some materials, usually presented as constant life¬ 
time diagrams , or master diagrams , such as those shown in Figure 5.39. If a designer is 
fortunate enough to find such data for a proposed material that match the operating condi¬ 
tions for the application of interest, of course these data should be used. If pertinent data are 
not available, an estimate of the influence of nonzero-mean stress may be made by utiliz¬ 
ing a mathematical model that approximates data of the type shown in Figure 5.38. 
Numerous successful models have been developed for this purpose. 32 Well known amongst 
these models and accepted as representing the experimental data are those due to Goodman 
(now called the modified Goodman relationship) and Gerber. The modified Goodman rela¬ 
tionship provides a simple linear estimate, whereas Gerber provides a parabolic estimate. 
Modified Goodman lines used to approximate the data are constructed on a plot of a a 


Figure 5.38 

Simulated high-cycle fatigue failure data 
showing the influence of mean stress. 
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32 See, for example, ref. 25, 26. 








Maximum stress, ksi Maximum stress, ksi Maximum stress, ksi 


Fluctuating Loads, Cumulative Damage, and Fatigue Life 259 


4.0 2.33 1.5 A = 1 0.67 0.43 0.25 0.11 0 


A = oo 

R = -1.0 
60 

50 
40 
30 
20 
10 
0 

4.0 2.33 1.5 A = 1 0.67 0.43 0.25 0.11 0 


A = oo 

R = -1.0 
150 

125 
100 
75 
50 
25 
0 

4.0 2.33 1.5 A = 1 0.67 0.43 0.25 0.11 0 


A = oo 

R = -1.0 
120 

100 

80 

60 

40 

20 

0 




AISI4340 

S n = 158 ksi, 

J? 

S y = 147 ksi 

RT, 2000 cpm 


- Unnotched 


- Notched 

K t = 3.3, p = 0.010 

J_I_I_I_L 

-120-100-80 -60 -40 -20 0 20 40 60 80 100 120 140 160 180 

Minimum stress, ksi 


-0.4 


- 0.2 


R = 0 


Ti-6 A1-4V 



Sol'n-treated 
and aged 
S ut =172 ksi, 
TYS = 158 ksi 
RT, 1850 cpm 

- Unnotched 

- Notched 

K t = 2.8 

j_i_2i_L 


-150-125-100-75 -50 -25 0 25 50 75 100 125 150 175 

Minimum stress, ksi 



-60 


7075-T6 
Wrought 
aluminum alloy 
S ut = 82 ksi, 
RT, 2000 cpm 

- Unnotched 

— Notched 
*,= 3.4 

J_I_I_L 

-50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 

Minimum stress, ksi 


Figure 5.59 

Master diagrams for alloys of 
steel, aluminum, and titanium. 
(From ref. 27, pp. 317, 322). 
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Figure 5.40 

Modified Goodman and Gerber diagram 
representing fatigue failure at a life of N 
cycles. 



versus a m by connecting the points S N on the a a axis to S u on the cr m axis in the tensile mean- 
stress region and passing a second line horizontally through S N in the compression 
mean-stress region. The Gerber parabola is constructed on a plot of a a versus cr m by con¬ 
necting the points S N on the cr a axis to S u , and ~S U on the a m axis in the tension-compression 
mean-stress region. These lines, shown in Figure 5.40, approximate the data illustrated in 
Figure 5.38 for a specified constant failure lifetime, N. Other failure lines would give simi¬ 
lar curves but would change since the magnitude of S N depends on failure life. Also it must 
be recognized that if \cr max \ or | cr min | exceed yield strength values of the material, failure by 
yielding will occur. Thus yielding failure lines may be plotted from a a = S yp on the cr a axis 
to a m = ±S yp on the a m axis, as shown in Figure 5.40 by dashed lines. 

Any combination of a a and a m that plots above either the solid or dashed lines in 
Figure 5.40 represents failure in a life shorter than N cycles, either by fatigue or by 
yielding. Writing an equation for line S N ~S U in the point-intercept form for the tensile 
mean-stress region, 

y = mx + b (5-61) 

which becomes 


Recalling that 



+ S N 


& a max ^ m 


and substituting into (5-62), 


C max 



— S N 


(5-62) 


(5-63) 


(5-64) 


or 


^mwc 



(5-65) 
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Defining the terms 


and 



R t 


(T 


m 


C max 


equation (5-65) may be rewritten as 


®max 


1 — m t R t 


(5-66) 

(5-67) 


(5-68) 


To clarify, a max in this equation represents the cyclic maximum stress that, with the pres¬ 
ence of mean stress a m , will cause failure in the same number of cycles as S N at zero-mean 
stress. To recognize that this is therefore a “strength parameter,” we set 


& max 


m S m 


(5-69) 


Further, cr max may not exceed S yp either, or yielding will occur. Summarizing then, for the 
tensile mean-stress region, 


5 'max—N 


1 — m t R t 


for CT m — 0 and Smax—N — S yp 


or 


where 


^ max—N 


S v , 


and 


for (r m > 0 


and 


c > c 

°max—N — °yt 


(5-70) 



R,= 


&m 

C max 


By a similar reasoning process for the compressive mean-stress region, a min represents 
the cyclic minimum stress that, with the presence of mean stress a m , will cause failure in 
the same number of cycles as S N at zero-mean stress. Again, \cr min \ may not exceed Is^l or 
yielding will occur. Thus for the compressive mean-stress region 


or 


where 


~N 


1 — m c R c 


for a m < 0 and S min - N > ~S yp 


Smin-N = -Syp for (J ,n ^ 0 and S min - N < ~S yp (5-71) 


and 


m r = 1 


Rr = 


Referring to Note 1 of Figure 5-40, and solving (5-62) from a different perspective, mod¬ 
ified Goodman expressions may be developed for an equivalent completely reversed cyclic 
stress , (J eq - CR . This stress cr eq _ CR is a calculated completely reversed cyclic stress that 
would produce failure in exactly the same number of cycles N as would the actual 
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nonzero-mean cyclic stress that has amplitude a a and mean a m . Thus any nonzero-mean 
cyclic stress may be transformed into an equivalent completely reversed cyclic stress of 
amplitude cr eq - CR that induces fatigue damage at the same rate cycle-by-cycle. Using this 
approach for the tensile mean-stress region 

<j ci 

®eq—CR h)r cr m — 0 and er max — Syp 

or 

®eq—CR Syp f^r — d and & max — Syp (5-72) 

and for the compressive mean-stress region 

&eq—CR ^ a ^m ~ d and l^mml ~ I Syp \ 

or 

V eq—CR I Syp | for (J m — 0 and \o"min\ — I Syp I (5-73) 

For the Gerber parabola, any combination of cr a and cr b that plots above either the 
dashed line or the solid curved line in Figure 5.40 again represents failure at a shorter life 
than N cycles, either by fatigue or yielding. Writing an equation for the Gerber parabola 
yields the form 

y = —ax 2 + b 

which becomes 

a a = -yyvl + S N (5-74) 

S u 

Again, referring to Note 1 of Figure 5.40 and solving (5-74) for an equivalent completely 
reversed cyclic stress, cr eq - CR , we have only at the tensile mean-stress region 

o\, 

a eq-CR = 

or 

® eq—CR Syp for rr m — 0 and cr max — Syp (5-75) 

The Gerber parabola is limited to tensile mean stresses because it incorrectly predicts a 
harmful effect for compressive mean stresses. 33 For compressive mean stresses (cr m < 0), 
equations (5-73) should be used. 

In all equations from (5-70) through (5-75), it should be noted that if the design ob¬ 
jective is infinite life, S f should be substituted for S N everywhere. 


>-'T 


for a n 


0 and a max < S yp 


Example 5.14 Design to Account for Nonzero Mean Stress 

A low-alloy steel link is to be made from a solid cylindrical bar subjected to an axial cyclic 
force that ranges from a maximum of 270 kN tension to a minimum of 180 kN compres¬ 
sion. The static material properties are S u = 690 MPa, S yp = 524 MPa, and elongation in 
50 mm of 25%. Calculate the diameter the link should be to provide infinite life. 


33 See ref. 28, page 430. 
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Solution 

Since the link is cyclically loaded and infinite life is desired, the fatigue limit Sf is the ma¬ 
terial property of primary interest; yielding should also be checked. Using the estimation 
guidelines, since S u < 200 ksi( 1379 MPa), 

S f = 0.5 (S u ) = 0.5(690) = 345 MPa 

Following the procedures of Example 5.12, and utilizing Table 5.3, the fatigue limit for the 
part may be estimated as 

Sf = kooSf = &oo(345) 

Since little information is given, the influencing factors of Table 5.3 may be estimated as: 

k gr =1.0 (from Table 5.3) 
k we =1.0 (no welding anticipated) 
kf = 1.0 (uniform shape specified) 
k sr = 0.65 (lathe-turned) 
k sz = 0.9 (from Table 5.3) 

k rs =1.0 (no information available; later review essential) 
kf r = 1.0 (no fretting anticipated) 

k cr =1.0 (no information available; later review essential) 

k sp =1.0 (from Table 5.3) 

k r = 0.9 (from Table 5.4 For R = 90) 


koo = (1.0)(1.0)(1.0)(0.65)(0.9)(1.0)(1.0)(1.0)(1.0)(0.9) = 0.5265 = 0.53 
S f = 0.53(345) = 182.9 MPa 

From the loading specification it may be noted that we have a case of nonzero-mean cyclic 
loading. Since stress is proportional to loading, 

Pm \{P max +P^) ~((270) + (—180)) 45 

°' m ~ ~A ~~ A _ A — ~A 


where A is the unknown cross-sectional area of the link. 

Since the mean stress is tensile, (5-67) may be used with 

S N =oo = Sf= 182.9 MPa 


690 - 182.9 

m t = ---= 0.735 

' 690 

a m PJA P„ 


so that (5-67) gives 


Since 


182.9 


45 

270 


= 0.167 


1 - m t R t 1 - (0.735)(0.167) 


= 208.5 MPa 


Smax-oo — 208.5 < Syp — 524 


( 1 ) 


Yielding does not occur and (1) is valid. 

To obtain the required diameter, we set the maximum applied stress cr max equal to the 
infinite-life strength for the specified nonzero-mean loading, giving 

max S max — oo 
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Example 5.14 
Continues 


or 


Solving for the diameter 


p 

1 max 

A 


4 p 


77 d 2 


208.5 


d = 


4P 


77(208.5) 


4(270 X 10 J ) 
77(208.5 X 10 6 ) 


= 0.0406 m = 0.041 mm 


Thus the required link diameter is 41 mm to provide infinite life. It should be noted, how¬ 
ever, that in practice a safety factor would also be imposed, and the required “safe” diam¬ 
eter would be somewhat larger than the 41 mm just calculated. 


Example 5.15 Design for Selected Reliability 

A one-inch-square bar is to be subjected to an axial cyclic force that ranges from a maxi¬ 
mum of 36,000 lb (tension) to a minimum of —22,000 lb (compression). The material is 
the wrought-steel alloy of Example 5.10 and the operating conditions are the same as those 
specified in Example 5.10. A fatigue strength reliability level of 99.9 percent is desired. 
How many cycles of operation would be expected before failure occurs? 

Solution 

Since the material and the operating conditions are exactly the same as for Example 5.10, 
and a fatigue strength reliability level of 99.9 percent is specified, the R = 99.9 percent 
curve of Figure E5.10 in Example 5.10 represents the appropriate estimate of the S-N curve 
for the one-inch-square bar. 

From the known loading and geometry, 


P mav 36,000 

2 = 36,000 psi 


1 max 

A 


( 1 ) 


Pn 

A 


mm - 22,000 

- — » -22,000 psi 


(l ) 2 


cr 


max + ^ 6,000 + (- 22 , 000 ) 


2 2 

^ 36,000 - (-22,000) 


7000 psi 
29,000 psi 


Since the mean stress is tensile, from (5-72) 


a a _ 29,000 

t ( 7000 \ 
S u V 76,000/ 


31,940 psi 


Also, 


d max = 36,000 < S yp = 42,000 


so yielding does not occur and cr eq - CR is valid. 

Since the mean stress is positive and yielding does not occur, Gerber’s equation could 
also be applied. Hence, using (5-75) 


a eq-CR ~ 


a 


a 


1 - 



29,000 
/ 7,000 \ 
V 76,000/ 


2 


29,248 psi 
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Figure E5.15 

Estimated S-N curve for fatigue 
strength reliability levels of 
50 percent (mean) and 
99.9 percent. 


Cycles to failure, N 


Referring to the S-N curve of Figure E5.15 for R = 99.9 percent, repeated below as 
Figure E5.15, and reading from the curve on the stress axis at 

s = o-eq-CR = 31,940 psi 

and 


S = cr eq - C R = 29,248 psi 

Then across to intersect the S-N curve, the intersections occur at N = 378,000 cycles, using 
Goodman and N = 827,000 cycles using the Gerber approach. 

The number of cycles to failure could also be obtained from the following equation 
since the S-N curve is taken as a straight line 


S = c + Mog(TV) 

Applying the conditions that at TV = 1, S = S u and at N = 10 6 , S'f = 0.5 S u , gives 


c = S,„ b = 


Sf - Su 


and 


logiN) 

N = K)(s-c)/fc 


From above, for R = 99.9 percent, S'f = 28,600 psi (see example 5.10) and, therefore, 

28,600 - 76,000 

c = 76,000 b = ---= -7900 

log(10 6 ) 







































































































266 Chapter 5 / Failure Theories 


Example 5.15 
Continues 


For S = 31,940 psi, 


N = 1 q (3 1940 -76000)/( - 790 °) = 378,000 cycles 


and for S = 29,250 psi 

N = io( 2925 °- 76000 V(-v 900) _ 8 27,000 cycles 

The appropriate fatigue failure curve to be used depends on how well the material property 
is known. 


Cumulative Damage Concepts and Cycle Counting 

Nearly all machines operate over a spectrum of speeds or loads, giving rise to a spectrum 
of alternating stress amplitudes and mean stresses, as illustrated for example in Figures 
5.25 and 5.26. Such variations in stress amplitudes and mean stresses make the direct use 
of standard S-N curves inapplicable because these curves are developed and presented for 
constant stress amplitude operation (and usually for zero-mean stress). Therefore it be¬ 
comes important to a designer to have available a theory or hypothesis, verified by exper¬ 
imental observations, that will permit good design estimates to be made for operation 
under conditions of spectrum loading using standard constant-amplitude S-N curves. 

The basic postulate adopted by all fatigue investigators is that operation at any given 
cyclic stress amplitude will produce fatigue damage , the seriousness of which will be related 
to the number of cycles of operation at that stress amplitude and also related to the total num¬ 
ber of cycles that would be required to produce failure of an undamaged specimen or part at 
that stress amplitude. It is further postulated that the damage incurred is permanent, and op¬ 
eration at several different stress amplitudes in sequence will result in an accumulation of to¬ 
tal damage equal to the sum of damage increments accrued at all the different stress levels. 
When the total accumulated damage reaches a critical value, fatigue failure occurs. 

Of the many damage theories proposed, 34 the most widely used is a linear theory first 
proposed by Palmgren in 1924 and later developed by Miner in 1945. This linear theory is 
referred to as the Palmgren-Miner hypothesis or the linear damage rule. The theory, based 
on the postulate that there is a linear relationship between fatigue damage and cycle ratio 
n/N , may be described using the S-N plot shown in Figure 5.41. 


Figure 5.41 

Illustration of spectrum loading where n t 



cycles of operation are accrued at each of &i 

the different corresponding stress levels 'f 
iTi , and the N- t are cycles to failure at ^ a 2 
each 07 . | 


Mean stress = 0 


< 


/\ / \ \ w \ 

n 2 n x N x n i _ 1 N 2 N i _ l u f N t 


34 See, for example, ref. 1, Ch. 8. 
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By definition of the S-N curve, operation at a constant stress amplitude will pro¬ 
duce complete damage, or failure, in Afi cycles. Operation at stress amplitude a 1 for a num¬ 
ber of cycles n u smaller than Afi, will produce a smaller fraction of damage, say D x . is 
usually termed the damage fraction. Operation over a spectrum of i different stress levels 
results in a damage fraction D f for each of the different stress levels in the spectrum. 
When these damage fractions sum to unity, failure is predicted; that is, failure is predicted 
to occur if (FIPTOI) 

D l + D 2 + • • • + A > 1 (5-76) 


The Palmgren-Miner hypothesis asserts that the damage fraction D, at any stress level 
cT t is linearly proportional to the ratio of number of cycles of operation n - L to the total num¬ 
ber of cycles N t that would produce failure at that stress level; that is, 


«i 




Nt 


(5-77) 


By the Palmgren-Miner hypothesis then, (5-76) becomes FIPTOI 

n i no H; 

—- + — + ••• + — > 1 

TVi N 2 N t 


(5-78) 


or FIPTOI 


n Hj 

s 1 (5-79) 

7=1 

This is a complete statement of the Palmgren-Miner hypothesis, or the linear damage 
rule. It has one sterling virtue, namely, simplicity; and for this reason it is widely used. 
Further, other much more complicated cumulative damage theories typically do not yield 
a significant improvement in failure prediction reliability. Perhaps the most significant 
shortcomings of the linear damage rule are that no influence of the order of application of 
various stress levels is recognized, and damage is assumed to accumulate at the same rate 
at a given stress level without regard to past history. Nevertheless, the Palmgren-Miner hy¬ 
pothesis may be used with good success most of the time. 


Example 5.16 Zero-Mean Spectrum Loading 

An alloy steel strut for an experimental aircraft application is fabricated from a supply 
of material that has an ultimate strength of 135,000 psi, yield strength of 120,000 psi, 
elongation of 20 percent in 2 inches, and fatigue properties under test conditions that 
match the actual operating conditions, as shown in the table of experimental results in 
Table E5.16. 

The cross-sectional area of the strut is 0.10 in 2 and buckling has been found not to be 
a problem due to the selected cross-sectional shape. In service, the strut is to be subjected 
to the following spectrum of completely reversed axial loads during each duty cycle: 

P a = 11,000 lb for 1000 cycles 
P b = 8300 lb for 4000 cycles 
P c = 6500 lb for 500,000 cycles 

This duty cycle is to be repeated three times during the life of the strut. 

Assuming fatigue to be the governing failure mode, would this strut be expected to 
survive all three duty cycles, or would it fail prematurely? 
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Example 5.16 
Continues 


110,000 

6,600 

105,000 

9,500 

100,000 

13,500 

95,000 

19,200 

90,000 

27,500 

85,000 

39,000 

80,000 

55,000 

75,000 

87,000 

73,000 

116,000 

71,000 

170,000 

70,000 

220,000 

69,000 

315,000 

68,500 

400,000 

68,000 

oo 


TABLE E5.16 Fatigue Test Data for Strut 
Material 

Stress Amplitude (psi) Cycles to Failure, N 


Solution 

Based on the data from Table E5.16, an S-N curve may be plotted as shown in Figure 
E5.16. 

Since the completely reversed axial loads are given, and the cross-sectional area of the 
strut is given as 0.10 in 2 , the completely reversed stress spectrum for each duty cycle may 
be calculated as 

P a 11,000 

(T a = — =-= 110,000 psi; n a = 1000 cycles 

a A q.10 F J 

P b 8300 

d h = — =-= 83,000 psi; n b = 4000 cycles 

/i U.1u 

P c 6500 

cr c ~ — =-= 65,000 psi; n c = 500,000 cycles 

c A 0.10 F c J 


Figure E5.16 

S-N curve for alloy steel data from Table 
E5.16. 
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Reading from the S-N curve of Figure E5.16 for completely reversed stresses a a , cr b , 
and cr c , respectively, 


N a = 6600 cycles to failure 

N b = 48,000 cycles to failure 

N c = oo (since cr c < Sf = 68,000) 


Furthermore, three such duty cycles must be sustained without failure. 

To determine whether failure would be predicted, the Palmgren-Miner linear damage 
rule formulated in (5-79) may be used. Hence FIPTOI 


or, FIPTOI 


it* i 

V ->1 

^ AT 


i=i N i 


_a_ 

Na N b N c _ 


Using the results above, FIPTOI 


1000 4000 500,000 

6600 + 48,000 + 


or FIPTOI 


3 [0.152 + 0.083 + 0] = 0.706 > 1 


Since this failure expression is not satisfied, the strut would be predicted to survive the 
three duty cycles. 


For the simple case of a stress spectrum composed of a sequence of uniaxial com¬ 
pletely reversed stresses of various amplitudes, the estimation of cumulative damage and 
prediction of failure are relatively straightforward, as illustrated by Example 5.16. If the 
stress spectrum is more complicated, as shown, for example, in Figures 5.25 or 5.26, the 
task of evaluating cumulative damage and even the task of counting the cycles in the spec¬ 
trum become much more difficult. Although numerous cycle counting methods have been 
devised, 35 the rain flow cycle counting method is probably more widely used than any 
other method and will be the only method presented in this text. 

To use the rain flow method, the stress-time history (stress-time spectrum) is plotted 
to scale on a stress-time coordinate system so that the time axis is vertically downward, 
and the lines connecting the stress peaks are imagined to be a series of sloping roofs, as 
shown, for example, in Figure 5.42. Several rules are imposed on “raindrops” flowing 
down these sloping roofs so that the rain flow may be used to define cycles and half-cycles 
of fluctuating stress in the spectrum. 

Rain flow is initiated by placing raindrops successively at the inside of each peak 
(maximum) or valley (minimum). 

The rules are as follows: 

1. The rain is allowed to flow on the roof and drip down to the next slope except that, 
if it initiates at a valley, it must be terminated when it comes opposite a valley more nega¬ 
tive than the valley from which it initiated. For example, in Figure 5.42 the flow begins at 


35 See, for example, ref. 1, p. 289. 
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Figure 5.42 

Example of rain flow cycle counting method. 
(After ref. 29.) 


Stress —► 



valley 1 and stops opposite valley 9, valley 9 being more negative than valley 1. A half¬ 
cycle is thus counted between valleys 1 and peak 8. 

2. Similarly, if the rain flow initiates at a peak, it must he terminated when it comes 
opposite a peak more positive than the peak from which it initiated. For example, in 
Figure 5.42 the flow begins at peak 2 and stops opposite peak 4, thus counting a half¬ 
cycle between peak 2 and valley 3. 

3. A rain flow must also stop if it meets the rain from a roof above. For example, in 
Figure 2.30, the flow beginning at valley 3 ends beneath peak 2. Using these rules, every 
part of the stress-time history is counted once and only once. 

4. If cycles are to he counted over the duration of a duty cycle or a “mission profile ” 
block that is to he repeated block after block, the cycle counting should be started by ini¬ 
tiating the first raindrop either at the most negative valley or the most positive peak, and 
continuing until all cycles in a complete block have been counted in sequence. This pro¬ 
cedure assures that a complete stress cycle will be counted between the most positive peak 
and most negative valley in the block. 

The nonzero-mean stress cycles may be converted to equivalent completely reversed cycles 
by utilizing (5-72) or (5-73). 

To determine the fatigue damage in each cycle associated with the equivalent completely 
reversed spectrum, an S-N curve for the material must be available or estimated using the 
techniques presented in Example 5.10 and the preceding discussion. 

Finally, damage is summed by utilizing the Palmgren-Miner linear damage rule 
(5-79) as illustrated in Example 5.16. 

It should be clear that for most real-life design cases the prediction techniques just de¬ 
scribed are practical only with the help of a digital computer program designed to carry 
forth the tedious cycle-by-cycle analyses involved. Many such programs are available. 36 

36 See, for example, ref. 15, App. 5A. 
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Example 5.17 Cycle Counting, Nonzero-Mean Loading, and Fatigue Life 

The stress-time pattern shown in Figure E5.17A is to be repeated in blocks. Using the rain 
flow cycle counting method and the S-N curve of Figure E5.17B, estimate the time in hours 
of testing required to produce failure. 

Solution 

The stress-time pattern of Figure E5.17A is first rotated so that the time axis is vertically 
downward, as shown in Figure E5.17C. Following the rules for rain flow cycle counting, 
the count is started at a minimum valley, as shown for raindrop (1) of the time-shifted 
block in Figure E5.17C. 

Data for each raindrop, read from Figure E5.17C, are recorded in Table E5.17. Values 
for cr eq - CR are calculated using (5-72) and (5-73). Values of N corresponding to each rain¬ 
drop’s cr eq - CR value are read from the completely reversed S-N curve of Figure E5.17B. 




Cycles to failure, N 


Figure £ 5 . J7A Figure E5.17B 

Stress-time pattern applied. S-N curve. 


Stress, ksi 


-50 -40 -30 -20 -10 0 10 20 30 40 



Figure E5.17C 

Stress-time pattern oriented for rain flow 
cycle counting. 


Time, minutes 
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Example 5.17 
Continues 


TABLE E5.17 Rain Flow Data for One "Block" of the Stress-Time Pattern 


Raindrop 

Number 

n , 

cycles 

^maxi 

ksi 

& min’ 

ksi 

^ nv 

ksi 

ksi 

a eq-CR’ 

ksi 

N, 

cycles 

1,4 @ \ cycle each 

1 

40 

-50 

-5 

45 

45 

6 X 10 3 

2,3 @ \ cycle each 

1 

20 

-10 

5 

15 

16.3 

OO 

5,8 @ \ cycle each 

1 

35 

10 

22.5 

12.5 

19.6 

OO 

6,7 @ \ cycle each 

1 

30 

20 

25 

5 

8.4 

OO 

9,10 @ \ cycle each 

1 

30 

0 

15 

15 

19.8 

OO 


Based on the Palmgren-Miner linear damage rule of (5-79), FIPTOI 



using Bf to denote the number of block to failure. 

From Table E5.17 it may be noted that only raindrops 1 and 4 will produce a nonzero 
cycle ratio. 

Thus FIPTOI 



> 1 


Using the equal sign, which corresponds to incipient failure, and reading numerical values 
from Table E5.17, 



1 


6 X 10 J 


1 


Bf = 6 X 10 3 blocks to failure 

At 1 block/minute, as indicated in Figures E5.17A and E5.17C, the predicted time to 
failure in hours, Hf, would be 

^6 X 10 3 blockslmin\/ 1 hr \ 


H f = 


"A 


block / V 60 min 


V 


= 100 hr to failure 


Multiaxial Cyclic Stresses 

Uniaxial cyclic stressing has been implicit in all of the discussion thus far in 5.6. Most real 
design situations, including rotating shafts, pressure vessels, power screws, springs, gears, 
flywheels, and many other machine elements, may involve multiaxial states of cyclic stress. 
Because of the complexities and costs in producing multiaxial fatigue failure data, only a 
limited body of such data exists. Consequently, no consensus has yet been reached on the 
best approach to prediction of failure under multiaxial fatigue stresses; however, various 
proposals have been discussed in recent books on fatigue analysis. 37 The approach adopted 
in this text involves the concept of equivalent stress to define a uniaxial equivalent to the 

37 See, for example, refs. 1, 19, 21, and 25. 
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actual multiaxial state of cyclic stresses, including both the alternating stress amplitudes and 
mean stresses. 


Fracture Mechanics (F-M) Approach to Fatigue 

In the introductory paragraphs of section 5.6 it was noted that fatigue is a progressive 
failure process that involves three phases: the initiation of crack nuclei, propagation of 
cracks until one reaches an unstable size, and finally, a sudden catastrophic separation 
of the affected part into two or more pieces. In the discussion of the S-N approach to 
fatigue design and analysis just completed, the three phases are lumped together in an 
S-N curve that represents a final failure locus for a part or specimen subjected to known 
fluctuating loads. The transitions from initiation to propagation to final fracture cannot 
be identified from an S-N curve, and the S-N approach does not require such informa¬ 
tion. Physically, however, the three phases have been well documented and the separate 
modeling of each phase in sequence has been under intense development. The results 
of these efforts have led to an approach to fatigue called the fracture mechanics (F-M) 
approach. An extensive database has now been established to support the F-M method 
and it has in recent years been very successfully used in the analysis and design of 
machine parts subjected to cyclic loads. A brief discussion of the F-M approach is in¬ 
cluded here, but much more complete discussions may be found in the recent literature 
on fatigue. 38 


Crack Initiation Phase 


The most widely accepted approach to the prediction of crack initiation life is the local 
stress-strain approach. While the details of the local stress-strain approach are compli¬ 
cated, the concepts are not. The basic premise is that the local fatigue response of the small 
critical zone of material at the crack initiation site, usually at the root of a geometrical dis¬ 
continuity, is analogous to the fatigue response of a small smooth laboratory specimen sub¬ 
jected to the same cyclic strains and stresses as the critical zone. This concept is illustrated 
in Figure 5.43. The number of cycles required to initiate a crack in the critical zone, N h is 
postulated to be equal to the number of cycles to produce failure of the small smooth spec¬ 
imen in a laboratory test under the same cyclic strains and stresses. Digital computer sim¬ 
ulation of the smooth specimen failure process then allows the prediction of N t if 
appropriate cyclic material response data are available. 

To proceed with the analysis, the prediction model must contain the ability to 

1. Compute local stresses and strains, including means and ranges, from the applied 
loads and geometry of the structure or machine part. 

2. Count cycles and associate mean and range values of stress and strain with each cycle. 

3. Convert nonzero-mean cycles to equivalent completely reversed cycles. 

4. Compute fatigue damage during each cycle from stress and/or strain amplitudes and 
cyclic materials properties. 

5. Compute damage cycle-by-cycle and sum the damage to give the desired prediction 
of N t . 

Many of these steps may be accomplished using methods already developed earlier in 
5.6. Rain flow cycle counting (Example 5.16) may be used for item (2), modified 
Goodman relationships (Example 5.15) for item (3), and Palmgren-Miner linear damage 
rule (Example 5.16) for item (5). Items (1) and (4) require further discussion. 

38 See, for example, refs. 1, 19, 21, and 25. 



Figure 5.43 

Smooth specimen analog 
of material at critical 
point in the structure. 
(See ref. 30.) 








274 Chapter 5 / Failure Theories 


To compute local stresses and strains from the external loading and geometry, a mod¬ 
ified version of the Neuber Rule 39 for estimating theoretical stress concentration factors 
leads to the expression 


AtxAe = 


(K f AS) 2 
E 


(5-80) 


where 


A a = local stress range 
Ae = local strain range 
A S = nominal 40 stress range 
Kf = fatigue stress concentration factor 41 
E = Young’s modulus of elasticity 


From other work, 42 an empirical expression for the cyclic stress-strain curve 43 satis¬ 
factory for most engineering metals has been developed as 


A e _ A a 
2 ~ 2 E 


A a 


W 


2k' 


(5-81) 


where k' and n' are material constants, the cyclic strength coefficient and cyclic strain¬ 
hardening exponent, respectively. Values for k' and n' may be found from the intercept and 
slope of an experimentally determined log-log plot of cyclic stress amplitude A<r/2 versus 
cyclic strain amplitude Ae/2. Values for many materials may be found in the literature 44 
With known values of material properties E, k', and n', the fatigue stress concentration fac¬ 
tor Kf, and readily calculable value of nominal stress range A S, equations (5-80) and 
(5-81) may be solved simultaneously to find values for local stress range A a and local 
strain range As. 

To compute the local fatigue damage associated with equivalent completely reversed 
stress and strain ranges for each cycle, it is necessary to have experimental failure data for 
strain amplitude versus cycles to failure (crack initiation), N h often plotted as strain am¬ 
plitude versus reversals to failure (crack initiation), 2 N h as illustrated in Figure 5.44. 

Noting from Figure 5.44 that total strain amplitude may be expressed as the sum of 
elastic strain amplitude plus plastic strain amplitude, each linear with cycles or reversals 
to failure (crack initiation) on a log-log plot, an empirical expression for total strain am¬ 
plitude versus cycles to crack initiation, N h has been developed as 45 

As (J 'f 

— = —{2N i ) b + S ' f {2N i y (5-82) 


The constants b and o-'f/E are the slope and one-reversal intercept of the elastic line in 
Figure 5.44, and the constants c and s' f are the slope and one-reversal intercept of the plas¬ 
tic line in Figure 5.44. The values of the cyclic material properties a' f (cyclic true fracture 
strength) and s' f (cyclic true fracture ductility) must be experimentally determined. With 


39 See, for example, ref. 1, p. 283. 

40 Stress calculated assuming no stress concentration effect; note that the use here of S for stress deviates from 
the usual policy in this text of reserving S to denote strength. 

41 See 5.3. 

42 See ref. 30, p. 7. 

43 Usually significantly different from the static stress-strain curve. See, for example, ref. 1, pp. 284-285. 
44 See, for example, ref. 31. 

45 See ref. 29. 
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Figure 5.44 

Schematic representation of elastic, plastic, 
and total strain amplitude versus fatigue life. 
(After ref. 29.) 


known cyclic material properties equation (5-82) may be used to calculate the number of 
cycles to crack initiation, N b for any given total strain amplitude, Ae/2. 

All five steps of the local stress-strain approach may therefore be completed if the 
cyclic material properties data are available and the applied nominal stress spectrum is 
known. 


Example 5.18 Fatigue Crack Initiation Life 

Experimentally determined values for the properties of a martensitic steel have been found 
to be S u = 1482 MPa, S yp = 1379 MPa, K Ic = 81 MPaVrn, e(50 mm) = 20%, 
k' = 1069 MPa, ri = 0.15, e' f = 0.48, <r' f = 2000 MPa, b = -0.091, and c = -0.60. A 
direct tension member made of this alloy has a single semicircular edge notch that results 
in a fatigue stress concentration factor 46 of 1.6. The rectangular net section of the member 
at the root of the notch is 9 mm thick by 36 mm wide. A completely reversed axial force 
of 70 kN amplitude is applied to the tension member. 

a. How many cycles would you estimate that it would take to initiate a fatigue crack at 
the notch root? 

b. What length would you estimate this crack to be at the time it is initiated by the cal¬ 
culation of part (a)? 

Solution 

a. The nominal stress amplitude S a may be calculated as 




70 000 


= 216 MPa 


A (0.009)(0.036) 

Hence the nominal stress range A S is given by 

AS = 2 S a = 2(216) = 432 MPa 

Using (5-75) 

[1.6(432 X 10 6 )] 2 

Ao-Ae =---= 2.31 MPa 

207 X 10 9 


( 1 ) 


46 See 5.3. 
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Example 5.18 
Continues 


Next, from (5-76), using the results from (1) 


A 8 
2 

A 8 
(As) 2 


A<t 


As 


+ 


2(207 X 10 9 )\A e) L2(1069 X 10 6 )VAsy. 


A <j 


As 


1/0.15 


2.3 X 10 b 
414 X 10 9 (As) 


+ 


2.3 X 10 b 


2138 X 10 6 (As) 


(1-1/0.15) 


§ 5.56 X 1(T 6 + 1.63 X 10 _20 (As) 

z 

Ae = Vl.ll X 10“ 5 + 3.26 X 10“ 20 (As)“ 5 ' 67 
This can be iterated to the solution 

As = 3.64 X 10 -3 m/m 


Then, from (5-77) 


3.64 X 10~ 3 _ 
2 

1.82 X 1(T 3 = 


2000 X 10 6 nnQ1 n , 

-—(27V/) _0-091 + 0.48(2^)“ 0-6 

207 X 10 9 

9.07 X 10“ 3 (A7 i )“ 0-091 + 0.317(7V;) -0-6 


This can be iterated to 


N t = 4.8 X 10 7 cycles to initiation 

b. There is no known method for calculating the length of a newly initiated fatigue crack. 
The length must either be measured from an experimental test or estimated from 
experience. Often, if no other information is available, a newly initiated crack is 
assumed to be approximately a t = 1.3 mm 


Crack Propagation and Final Fracture Phases 

As discussed in 5.5, the concepts of linear elastic fracture mechanics may be employed to 
predict the size of a crack in a given structure or machine part that will, under specified 
loadings, propagate spontaneously to final fracture. This critical crack size , a cn may be de¬ 
termined, for example, by solving (5-52) for crack size, which then (by definition) be¬ 
comes a cr when the applied stress a is given its maximum value. Thus, from (5-52) 


O'er _ 

77 


K r 


CCTyy 


(5-83) 


Further, if plane strain conditions are met , as defined by (5-53), K c becomes equal to 
the plane strain fracture toughness K Ic , and (5-83) becomes 


a 


cr 


K„ 


C(T n 


(5-84) 


A fatigue crack that has been initiated by cyclic loading, or any other preexisting flaw 
in the material, may be expected to grow under sustained cyclic loading until it reaches the 
critical size, a cr , from which it will propagate spontaneously to catastrophic failure in ac¬ 
cordance with the laws of fracture mechanics. Typically, the time for a fatigue-initiated 
crack to grow to critical size is a significant portion of the life of the machine part. Thus, 
not only is it necessary to understand the crack initiation phase and the definition of criti¬ 
cal crack size but an understanding of the growth of a crack from initial size a t to critical 
size a cr is also essential. 
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Figure 5.45 

Fatigue crack propagation behavior of various 
steels. (Reprinted from ref. 34 by permission of 
Pearson Education, Inc., Upper Saddle River, 
N.J.) 


The crack growth rate, da/dN, has been found to be related to the range of stress- 


intensity factor, A K, where 


A K = CA(t\/ttci 

(5-85) 

and the range of stress is given by 


O' max Omin 

(5-86) 


Most crack propagation data are presented as log-log plots of da/dN versus A K. Figure 
5.45 shows such a plot for various types of steels. The study of such data has led to many 
different empirical models 47 for the prediction of crack growth rate da/dN as a function of 
A K. One such model, developed by Paris and Erdogan, 48 has become widely known as the 
Paris Law , 


J 7 = C PE {XK) n (5-87) 

where n is the slope of the log-log plot of da/dN versus A K, as shown in Figure 5.45, and 
C PE is an empirical parameter that depends upon material properties, cyclic frequency, 
mean stress, and perhaps other secondary variables. 

47 See, for example, ref. 32. 48 See ref. 33. 
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If the parameters C PE and n are known for a particular application, the crack length, 
a N , resulting from the application of N cycles of loading after the crack is initiated, may be 
computed from the expression 

n p 

a N = a, + 2 C PE (SK *) (5-88) 

j= 1 

or 



where a t is the length of a newly initiated crack and N p the number of loading cycles in the 
propagation phase following crack initiation. For complicated spectrum loading histories, 
these computations require block-by-block or cycle-by-cycle analyses, making the use of 
modern digital computer systems mandatory for the solution of most practical design prob¬ 
lems. It should also be recognized that various other phenomena may take place in the zone 
of material around the tip of a growing crack to influence the crack growth prediction 
accuracy. 49 These may include crack growth retardation or acceleration due to plastic 
zones produced by preceding cyclic load histories, the influence of nonzero-mean stresses, 
and crack size effects. Typically, a designer would utilize (5-88) or (5-89) to obtain 
preliminary estimates of crack propagation life, but would enlist the help of a fracture 
mechanics specialist to improve prediction accuracy by including modifications to account 
for such factors as loading sequence, environment, frequency, multiaxial states of stress, 
and determination of applicable A K values in view of plasticity at the crack tip. 


Example 5.19 Fatigue Crack Propagation Life 

The martensitic alloy steel tension member of Example 5.17 is to be subjected to a contin¬ 
uation of completely reversed cyclic axial force of 70 kN amplitude, following the initia¬ 
tion of a fatigue crack. It may be assumed that the length of the newly initiated fatigue crack 
is 1.3 mm. Further, it may be assumed that the stress concentration effects are negligible for 
the propagating crack (because the initiated crack tip probably extends through most of the 
zone of stress concentration). How many cycles of continued loading would you estimate 
could be applied before catastrophic fracture would occur? 

Solution 

Since the applied axial force is completely reversed, the maximum and minimum normal 
stresses are 


Fa _ 70 000 

± A ~ ± (0.009)(0.036) 


±216 MPa 


Using (5-50) to check the plain strain conditions 

B = 0.009 > 2.5 (—) =2.5(±U) = 0.0086 
\Sy„J V 1379/ 

Since this expression is satisfied, plain strain conditions exist and the critical crack size 
may be determined from (5-79). Thus, 


1 

K Ic ' 

2 

1 

81 

2 0.0448 

7 T 

- Ccr max _ 


7 T 

_216C_ 

c 2 


The parameter C may be evaluated from Figure 5.19 using the methods of Example 5.9, 
but must be iterated together with (1) above since the crack length a cr must be known to 


49 See, for example, ref. 1, pp. 293-304. 
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find C, and C must be known to find a cr . To illustrate the iteration, assume, for example, 
a cr = 15 mm. Then from Figure 5.19, 

alb = 15/36 = 0.42 


and 

C(1 - 0.42) 3/2 = 0.97 C = 2.2 

From (1) 

0.0448 

a cr =-— = 0.0093 m = 9.3 mm 

(2.2) 2 

Iterating, try a cr = 12.5 mm, giving alb = 12.5/36 = 0.35 and C = 0.98/(l - 0.35) 3/2 = 
187. This results in 

0.0448 

a cr = -- = 0.0128 m = 12.8 mm 

(1.87) 2 

This is considered close enough to the assumed value of a cr = 12.5 mm. Corresponding to 
this critical crack length, we use C cr = 1.87. 

For the newly initiated crack of length a t =1.3 mm, the corresponding C, may be 
found from Figure 5.19, with ajb = 1.3/36 = 0.036 

C, = 1.07/(1 - 0.036) 3/2 = 1.13 

Since this is a martensitic steel, from Figure 5.45 the applicable empirical crack growth 
model must be converted into SI and is 

— = 3.03 X 1(T 10 (A K) 2 - 25 (2) 

dN 

where A K is in MPa Vm and da/dN is in meter per cycle. 

Next we establish A a using the maximum and minimum normal stresses 

Ao- = 216 - (-216) = 432 MPa 

Substituting (5-80) into (2) above 

= 3.03 X 10 _10 (CAcrV / T7a) 2 ' 25 


or 


da 
_ 1.125 


= 3.03 X 10“ 10 (CAuV7r) 2 - 25 JAf 


Integrating both sides 

r a r r =0.0125 


/ 


da 


a.=0.0013 a 


1.125 


« 3.03 X 10 


-10 


1.13 + 1.87 


(432)Vi7 


2.25 rN p 


dN 


(3) 


Note that C is taken as the average of C { and C cr . A more accurate solution could be found 
by partitioning the crack growth into smaller increments and separately integrating each 
increment, summing the results to obtain N p . Evaluating (3) 


—0.125 


0.125 


0.0125 


= 0.00233/V n 


0.0013 


N p = 1.94 X 10 3 


N = 1.94 X 10 3 propagation cycles to final fracture 
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Design Issues in Fatigue Life Prediction 

Using the F-M approach to estimate the total fatigue life of a proposed design configura¬ 
tion is simple in concept. The total fatigue life to failure N f is the sum of the initiation life 
N t plus the propagation life N p . That is, 

N f = N t + N p (5-90) 

The initiation life may be calculated using the local stress-strain approach, as illus¬ 
trated in Example 5.18. The propagation life may be calculated using an appropriate crack 
growth model together with LEFM estimates of critical crack size, as illustrated in 
Example 5.19. As a practical matter, however, several important additional issues must be 
addressed when attempting to use (5-90) as a design tool. These include the following: (1) 
determining or specifying the size of an initiated crack corresponding to N t cycles; (2) ac¬ 
counting for geometric stress concentration effects and stress gradients; (3) accounting for 
strength gradients, especially as associated with metallurgical or mechanical surface treat¬ 
ments; (4) accounting for residual stress fields; (5) accounting for multiaxial states of 
stress, three-dimensional effects, and others. Most of these issues remain research topics, 
requiring designers to consult with fracture mechanics specialists and/or make appropriate 
simplifying assumptions when utilizing (5-90). In the final analysis, it is essential to con¬ 
duct full-scale fatigue tests to provide acceptable reliability. 


Example 5.20 Estimating Total Fatigue Life 

Referring to the solutions of Examples 5.18 and 5.19, estimate the total fatigue life of the 
martenstic steel tensile member subjected to a completely reversed cyclic axial force of 70 
kN amplitude. 

Solution 

The total fatigue life may be estimated using (5-90). With the results from Examples 5.18 
and 5.19 

N f = N { + N p = 4.8 X 10 7 + 1.94 X 10 3 = 4.8 X 10 7 cycles 
For this member, the fatigue life Nf is dominated by the crack initiation phase. 


Fatigue Stress Concentration Factors and 
Notch Sensitivity Index 

Unlike the theoretical stress concentration factor K t , the fatigue stress concentration factor K f 
is & function of the material , as well as geometry and type of loading. To account for the in¬ 
fluence of material characteristics, a notch sensitivity index q has been defined to relate the 
actual effect of a notch on fatigue strength of a material to the effect that might be predicted 
solely on the basis of elastic theory. The definition of notch sensitivity index q is given by 


where Kj = 
K t = 
q = 


fatigue stress concentration factor 

theoretical stress concentration factor 

notch sensitivity index valid for high-cycle fatigue range 


The magnitude of q ranges from zero for no notch effect to unity for full notch effect. The 
notch sensitivity index is a function of both material and notch radius, as illustrated in 
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Notch radius r, in. 


Figure 5.46 

Curves of notch sensitivity index 
versus notch radius for a range of 
steels and an aluminum alloy 
subjected to axial, bending, and 
torsional loading. (After ref. 15; 
reprinted with permission of the 
McGraw-Hill Companies.) 


Figure 5.46 for a range of steels and an aluminum aloy. For finer-grained materials, such 
as quenched and tempered steels, q is usually close to unity. For coarser-grained materials, 
such as annealed or normalized aluminum alloys, q approaches unity if the notch radius 
exceeds about one-quarter inch. In view of these facts it is tempting to recommend the use 
of Kf = K t as a simplifying assumption. Doing so, however, would ignore several impor¬ 
tant notch sensitivity effects, including: 

1. Under fatigue loading, an alloy steel with superior static properties will often be found 
not to have superior fatigue properties when compared to a plain carbon steel, because 
of the difference in notch sensitivities. 

2. There is a tendency to improperly assess the effects of tiny scratches and cavities un¬ 
less notch sensitivity effects are recognized. 

3. Serious errors in applying the results from models to large structures may be made if 
notch sensitivity effects are not recognized. 

4. In critical design situations, inefficiencies may accrue if notch sensitivity effects are 
not considered. 

Based on (5-91), an expression for fatigue stress concentration factor may be written as 

K f — q(K t — 1) + 1 (5-92) 

The theoretical elastic stress concentration factor K t may be determined, on the basis of 
geometry and loading, from handbook charts such as those depicted in Figures 5.4 
through 5.12. The notch sensitivity index q may also be read from charts, such as the one 
shown in Figure 5.46. 

For uniaxial states of cyclic stress it is sometimes convenient to use K f as a “strength 
reduction factor” rather than as a “stress concentration factor.” This may be done by di¬ 
viding the fatigue limit by ^rather than multiplying the applied nominal cyclic stress times 
K f . Although conceptually it is more correct to think of K f as a stress concentration factor, 
computationally it is often simpler to use K f as a strength reduction factor when the cyclic 
stresses are uniaxial. For multiaxial states of stress, however, K f must be used as a stress 
concentration factor. 

The fatigue stress concentation factor (or strength reduction factor) determined from 
(5-92) is strictly applicable only in the high-cycle fatigue range (lives of 10 1 2 3 4 5 -10 6 cycles 
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Figure 5.47 

S-N curves for notched and unnotched 
specimens subjected to completely 
reversed axial loading. (After ref. 16.) 



Cycles to failure 


and greater). It has earlier been noted that for ductile materials and static loads, effects of 
stress concentration may usually be neglected. In the intermediate and low-cycle life range 
(from a quarter-cycle up to about 10 5 -10 6 cycles), the stress concentration factor increases 
from unity to K f , so the notched and unnotched S-N curves tend to converge toward the 
quarter-cycle point A, as shown in Figure 5.47. 

Estimates for fatigue stress concentration factor can be made by constructing a 
straight line from the ultimate strength plotted at a life of 1 cycle to the unnotched fatigue 
strength divided by K f plotted at a life of 10 6 cycles, on a semilogarithmic S-N plot. 50 The 
ratio of unnotched to notched fatigue strength at any selected intermediate life may be used 
as an estimate of fatigue stress concentration factor for that life. 

Finally, experimental investigations have indicated that for evaluating fatigue of duc¬ 
tile materials , the fatigue stress concentration factor for any nonzero mean cyclic state of 
stress should be applied only to the alternating component of stress (and not to the steady 
component). In evaluating the fatigue loading of brittle materials, the stress concentration 
factor should be applied to both the alternating and the steady component. 


Example 5.21 Fatigue Life Prediction Under Uniaxial Stress, Including Stress 
Concentration 

A 60-mm-wide by 2.5-mm-thick rectangular plate of annealed 1040 steel has a 6-mm- 
diameter hole drilled through, as shown in Figure E5.21. The properties of the 1040 steel 
are S u = 372 MPa, S yp = 330 MPa, e(50 mm) = 50%, and S f = 186 MPa. The plate is to 
be subjected to a completely reversed alternating direct force of 8 kN. Buckling is not a 
problem. What life cycle would you predict for the bar if the not notched S-N curve of 
Figure 5.47 is valid? 

Solution 

The actual stress amplitude at the critical point adjacent to the hole is 

_ ( P a \ _ ( 8000 

((Ta)act ~ K ^ a)nom ~ K \A net ) ~ ^0.06(0.0025) 




t = 2.5 mm / 



d = 6mm \ ) 

\ 

b = 60 mm 



j = 53.3/^MPa 


Figure E5.21 

Steel bar with stress concentration p a = 8 kN 
subjected to cyclic loading. ^— 


P = 8 kN 
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From (5-92) 


K f = q(K t - 1) + 1 

Using Figure 5.46 with S u = 372 MPa = 54 ksi and r = 3 mm ~ 0.125 in, 

q = 0.76 

From Figure 4.9 (b), with bid = 6/60 = 0.10, 


Therefore, 


K t = 2.7 


K f = 0.76(2.7 - 1) + 1 = 2.29 
The actual stress is, therefore, 


(Oacf = 53.3 K f = 53.3(2.29) = 122 MPa 
Calculating the ordinate ratio 


(^a)act 

S u 


-= 0.378 = 37.87% 

372 


And using it to read into the “not notched” S-N curve of Figure 5.47, the predicted cyclic 
life of the plate is infinite. 51 It should be recognized, however, that there is virtually no 
margin of safety. It would be prudent, therefore, to impose an appropriate safety factor or 
reliability assessment before going further. 


5.7 Multiaxial States of Cyclic Stress and Multiaxial 
Fatigue Failure Theories 


An extensive discussion of the prediction and prevention of fatigue failure and fatigue life 
under uniaxial states of cyclic stress was presented in 5.6. As mentioned in that discussion, 
most real design situations involve fluctuating loads that produce multiaxial states of cyclic 
stress. A consensus has not yet been reached on the best approach to prediction of failure 
under multiaxial states of cyclic stress. However, one technique used for ductile materials 
subjected to fluctuating multiaxial stresses is to utilize a cyclic adaptation of the equiva¬ 
lent stress expression of (5-45). For brittle materials subjected to fluctuating multiaxial 
stresses, the principal normal stress expressions are utilized. 

An equivalent uniaxial alternating stress expression <r eq - a for ductile materials, based 
on the equivalent stress equation of (5-47), becomes 


O eq—a yj r^_(o x — a Vy—a) (Vy—a O z — a ) "f" (o z — a &x—a) ] 3(t xy—a fyz~a ^xz~a ) 

(5-93) 

An equivalent uniaxial mean stress expression cr eq _ m for ductile materials, based on the 
equivalent stress equation of (5-47), is 


o'. 


eq—m 


2 ^ 


Oy—m) (Oy — m O z — m ) + (cr z — n 


O r 


„) 2 ] + 3(r 


xy—m 


+ t 2 + r 

1 * yz—m ' * xz 


—m) 


(5-94) 


5 'Because the ordinate ratio 0.42 lies below the S-N curve. 
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or, for plane stress, 

O'eq—a K — a 4" *Jy—a®x—aP'y—a 4" ^ 7 xy—a (^”^5) 

Ceq—m ^&x~m 4 " &y—m &x-mP'y—m 4 " xy—m (5-96) 

The equations (5-93) through (5-96) may be substituted into any of the expressions 
discussed in 5.6 as equivalents for their uniaxial counterparts, a a and cr m . For example, in 
using (5-72), if the loading conditions produce a multiaxial state of cyclic stress the equa¬ 
tions become 

&eq-a 

® eq-CR i / /n \ for CT e q_ m — 0 and & max — ^yp (^”97) 


where 


®eq-CR Syp 


for cr e q_ m — 0 and cr max ~ ^yp (^ 98) 


®max ®eq-a 4" ^eq-m (5-99) 

For brittle materials subjected to multiaxial states of cyclic stress, if the convention 
0 -j > <t 2 > <r 3 is adopted, the expressions a Xa and cr lm may be substituted into any of the 
uniaxial procedures discussed in 5.6. With brittle materials, for example, equations (5-73) 
would become, under multiaxial states of cyclic stress, 


a eq-CR ~ 


(J la 

1 (^"l m/^u) 


for a , m > 0 and cr max ^ S u (5-100) 


Example 5.22 Fatigue Failure Prediction Under Fluctuating Multiaxial Stresses 

A power transmission shaft of solid cylindrical shape is to be made of 1020 hot-rolled steel 
with S u = 65,000 psi, S yp = 43,000 psi, e = 36 percent elongation in 2 in, and fatigue prop¬ 
erties as shown for 1020 steel in Figure 5.31. The shaft is to transmit 85 horsepower at a 
rotational speed of n = 1800 rpm, with no fluctuations in torque or speed. At the critical 
location, midspan between bearings, the rotation shaft is also subjected to a pure bending 
moment of 1500 in-lb, fixed in a vertical plane by virtue of a system of symmetrical 
external forces on the shaft. If the shaft diameter is 1.0 in, what operating life would be 
predicted before fatigue failure occurs? 

Solution 

From the horsepower equation (4-39), the steady shaft torque is 


63,025 (hp) _ 63,025(85) 
n 1800 


2976 in-lb (steady) 


And from the problem statement, the bending moment is completely reversed (due to shaft 
rotation), giving 


M= 1500 in-lb (completely reversed) 

Since the maximum shearing stress due to torque T occurs at the surface, and cyclic 
bending stresses range from maximum to minimum and back to maximum at the surface 
with each rotation, all midspan surface points are equally critical. 

A typical volume element at a midspan critical point is shown in Figure E5.22. The 
steady torsional shearing stress may be calculated from (4-37) as 


_ Ta _ 16 T 
J 7Td 3 
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z Figure E5.22 

State of stress at a typical midspan critical 
point. 



and cyclic bending stress from (4-5) as 

Me 32 M 


Noting that T max = T min = T m = T = 2976 in-lb, M max = +1500 in-lb, M min = -1500 in¬ 
lb, and M m = 0, the alternating and mean components stresses are 

16(2976) 

^ xy—m — 3 — 15,157 psi, T xy - a — 0 psi 

TT{\y 

32(1500) 

Cx-m = Opsi, (T x - a = -— 1 — = 15,279 psi 

7r(iy 

Since e = 36 percent, the material is ductile and the equivalent stress concepts of 
(5-95) and (5-96) may be used to write 


a eq - a = V(15,279) 2 + 0 2 = 15,279 psi 
a eq - m = Vo 2 + 3(15,157) 2 = 26,253 psi 


Using the expressions in (5-99), the condition of cr max may be checked by calculating 
a max = 15,279 + 26,253 = 41,532 < 43,000 psi 
Thus, (5-97) is valid, and gives 


15,279 

^- CR = -26^53 = 25 ’ 631psi 

1 65,000 

Checking Figure 5.31, this value, cr eq . CR = 25,631 psi, corresponds to infinite life for 1020 
steel. No strength-influencing factors such as surface finish, corrosion, and so forth, have 
been considered. A more accurate prediction could be made by considering such factors, 
as was done in Example 5.13. 


Example 5.23 Fatigue Life Predication Under IVIultiaxial Stresses, Including Stress 
Concentration Effects 

Refer to Figure E5.23A, in which a torsionally oscillating shaft of 1.25-in diameter has 
a 0.25-inch diametral hole all the way through it. By the way the shaft is loaded, it is 
subjected to a released torsional moment of 8300 in-lb, and an in-phase released cyclic 
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Example 5.23 bending moment, in the plane of the through-hole axis, of 3700 in-lb. If the shaft is made 

Continues of 4340 steel with S u = 150,000 psi, S yp = 120,000 psi, e = 15 percent in 2 in, and the 

fatigue properties shown in Figure E5.23B, how many torsional oscillations would you 
expect could be completed before fatigue failure of the shaft takes place? (Assume that 
the critical point for bending and torsion coincide.) 



Figure E5.23A 

Oscillating shaft with through hole. 



Cycles to failure 


Figure E5.23B 

Uniaxial fatigue properties of 4340 steel used in the torsion- 
ally oscillating shaft. 


Figure E5.23C z 

State of stress at the critical point. 



Solution 

The critical points for bending and torsion are assumed to be at the same location, adja¬ 
cent to the edge of the hole at top and bottom surface elements. The state of stress at the 
critical point is shown in Figure E5.23C. 

From Figure 5.6, the nominal bending and torsional stresses as well as bending and 
torsional stress concentration factors may be found. The nominal bending stress is 

M M 

(T n0 m = -o-^ = -7-- = 7.9M 

7rP 3 _ drf 77(1. 25) 3 0.25(1.25) 2 

32 6 32 6 

and the nominal torsional stress is 


Tnom 


T 

7 tD 3 _ dD ^ 

^L6~ ~ 6 ~ 


T 

-r-y = 3.14T 

77(1. 25) 3 0.25(1.25) 2 

16 6~ 


025 
D 1.25 


Referring again to Figure 5.6, for 
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It may be found that 


(K t ) tor = 1.50 


and 


Also, from Figure 5.46 


and 


Then from (5-92) 


and 


Then 


and 


(K t ) b end = 2.03 


q tor = 0.94 


Qbend 0.92 


K f _ tor = 0.94(1.50 - 1) + 1 = 1.47 

K f - bend = 0.92(2.03 - 1) 4 1 = 1.95 


r xy = 1.47(3.147) = 4.62T 


a x = 1.95(7.90M) = 15.41AT 

Since the shaft is subjected to a released torsional moment of 8300 in-lb and a released 
cyclic bending moment of 3700 in-lb, the maximum and minimum values are 

M max = 3700 in-lb, M min = 0 
^max = 8300 in-lb, T min * 0 

The alternating and mean values are 

3700 

M a — M m — ~y~ = 1850 in-lb 
8300 

T a = T m = — = 4150 in-lb 

Based on the above equations, the alternating and mean component stresses are 

d x - a — cr x - m = 15.41(1850) = 28,509 psi 
Txy-a = Txy-m = 4.62(4150) = 19,173 psi 

Since e = 15 percent, the material is ductile, and the equivalent stress concepts of (5-95) 
and (5-96) may be used to give 

<r eq - a = V(28,509) 2 + 3(19,173) 2 = 43,767 psi 


And 

cr eq - m = V(28,509) 2 + 3(19,173) 2 = 43,767 psi 
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Example 5.23 Using these expressions in (5-99), the condition on cr max may be checked by calculating 

Continues 

o-max = 43,767 + 43,767 = 87,534 psi 
Based on the validity criteria 


And 


a max = 87,534 < S yp = 120,000 psi 
If (5-72) is chosen, the result is 

43,767 

<7e?_CR 1 - (43,767/150,000) 61 ’ 800 P S1 

If (5-75) is chosen, the result would be 


a eq~CR 


43,767 

-- = 47,840 psi 

1 - (43,767/150,000) 2 


Since the fatigue limit in Figure E5.23B is shown as Sf = 75,000 psi, both values calcu¬ 
lated for cr eq _ CR would be predicted to give infinite life (an infinite number of torsional os¬ 
cillations) for this part. Strength-influencing factors such as surface finish, corrosion, and 
so on, should be reviewed, however, to improve the life assessment. (See Example 5.13.) 


Example 5.24 Fatigue Life Improvement Resulting from Shot-Peening (Favorable 
Residual Stresses) 

A cylindrical support bar is to be subjected to a released cyclic tensile load of P max = 23,000 
lb. The critical section in the bar is at the root of a circumferential groove of semicircular 
profile, having the dimensions shown in Figure E5.24A. The material is to be a wrought 
carbon steel alloy with S u = 76,000 psi, S yp = 42,000 psi, and e (2 inches) =18 percent 
when in hot rolled condition, and S u = 70,000 psi, S yp = 54,000 psi, and e (2 inches) =15 
percent when in cold drawn condition. Fatigue properties are defined by the S-N curve 
from Example 5.10, and repeated here as Figure E5.24B. 

a. Using a fatigue strength reliability level of 99.9 percent, estimate the life of the bar 
under the conditions given, if the material is in hot-rolled condition. 

b. Estimate the improvement in life (if any) that might be expected from the hot-rolled 
bar if the groove is properly shot-peened, assuming that the shot-peening process 
results in a compressive residual surface stress of about 50 percent of the strain- 
hardened yield strength. 


Figure E5.24A 

Cyclically loaded cylindrical steel support 
bar with circumferential groove of semi¬ 
circular profile. 
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Figure E5.24B 

Estimated S-N curves for 
fatigue strength reliability 
levels of 50% and 99.9% 
(repeat of Figure E5.10A). 
Note that Sf for R = 50% 
is 38,000 psi and Sf for 
R = 99.9% is 28,600 psi. 


Solution 

a. For the released loading applied, with P max = 23,000 lb, the nominal cyclic stress 
ranges from 

( & min ) nom ^ 


to 


P max 23,000 

( \ _ max _ * 

a max)nom ~ A — , , 2 

A net 77(1.00) 


29,285 psi 


giving nominal values of stress amplitude and tensile mean stress of 

29,285 - 0 

( O- a ) n0 m = -^- = 14 ’ 642 P S1 


and 


29,285 + 0 

(o- m )nom = - 2 - = 14,642 psi 


Using methods developed in 5.3, since the material is ductile, 

(Va)act Kf(Va)nom 


and 


(,®~m)act ( ^ m ) nom 

Then from (5-91) 

K f — q(K t — l) + l 

where q may be read from Figure 5.11 with r = 0.06 and S u = 76,000 psi, as 

q = 0.74 

From Figure 5.5(b), with r/d = 0.06 and D/d = 1.12, K t may be read as 

K t = 2.18 
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Example 5.24 
Continues 


giving 


From this, 


and 


K f = 0.74(2.18 - 1) + 1 § 1.87 
(o- a ) act = 1.87(14,642) = 27,381 psi 


( o- m ) ac t = 14,642 psi 

Therefore, effective values for cr max and a min are 

(<r max ) ej f= 27,381 + 14,642 = 42,023 psi 


and 


(■ <r min )e ff = 14,642 - 27,381 = -12,739 psi 


Now, using (5-72), 


Veq-CR ~ 


(.Cajact 

1 C m 


27,381 
14,642 
~~ 76,000 


33,915 psi 


Taking this equivalent completely reversed stress amplitude to the 99.9 percent relia¬ 
bility S-N curve of Figure E4.20B, the estimated life is read as 

N = 2 X 10 5 cycles 

b. If the groove were shot-peened to give a residual compressive stress of 50 percent of the 
strain-hardened yield strength, the compressive residual stress at the surface would be 

a res — —0.5(54,000) = —27,000 psi 

Superposing the residual stress on the effective operational stresses at the notch 
root, we get 


{o’max) S p-eff= 42,023 + (-27,000) = 15,023 psi 
and 

{(T m in)sp-eff= -12,739 + (-27,000) = -39,739 psi 
so 

15,023 - (-39,739) 

i.(Ta)sp-eff = - 2 -= 27,381 psi 

and 

15,023 + (-39,739) 

^ — 12,358 psi 

For the shot-peened groove, then, since mean stress is compressive, using (5-73), 
Veq-CR = Wa)sp-eff = 27,381 psi 
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Taking the equivalent completely reversed stress amplitude to the 99.9 percent re¬ 
liability S-N curve of Figure E5.24B, the estimated life is read as 

N sp = oo cycles 


because 


O'eq-CR = 27,381 < S f (R = 99.9%) = 28,600 psi 

Shot-peening the groove in this case is estimated to increase fatigue life of the part 
from about 200,000 cycles to infinite life. Shot-peening is widely used to improve 
fatigue failure resistance. 



